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2 Abstract. The scaling of plate-tectonic convection is investigated by simulating thermal 

3 convection with pseudoplastic rheology and strongly temperature-dependent viscosity. The 

4 effect of mantle melting is also explored with additional depth-dependent viscosity. Heat-flow 

5 scaling can be constructed with only two parameters, the internal Rayleigh number and 

e the lithospheric viscosity contrast, the latter of which is determined entirely by rheological 

7 properties. The critical viscosity contrast for the transition between plate-tectonic and 

s stagnant-lid convection is found to be proportional to the square root of the internal Rayleigh 

9 number. The relation between mantle temperature and surface heat flux on Earth is discussed 

10 on the basis of these scaling laws, and the inverse relationship between them, as previously 

11 suggested from the consideration of global energy balance, is confirmed by this fully dynamic 

12 approach. In the presence of surface water to reduce the effective friction coefficient, the 

13 operation of plate tectonics is suggested to be plausible throughout the Earth history. 
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1. Introduction 



15 



Simulating mantle convection wi th plate tectonics in a fully dynamic manner has become 



popular in the last decade or so [e.g., 



Bercovici 



2003I1 . and quite a few studies have been 



17 published addressing a variety of problems, including t he significance of 3-D spheric al 



s geometry wichards et al. 



2001 



van Heck and Tackley, 



9 the ro le of history-dependent rheol ogy [e.g., 



20 



2008], the initiation of subductio n jSolomatov 



21 to other terrestrial planets [e .g.. 



Landuyt and Bercovici . 



Tackley, 



2004 



Lenardic et al. 



2008 



200ol: 



Foley and Becker, 



Ogawa, 



Gurnis et al. 



2004; 



2003 
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2009], 



Landuyt et al. . 



2004] , and appli cations 



O'Neill and Lenardid 



2007; 



2009] . Around the same time, interests in the i nitiation and evolutio n 
23 of plate tecton i cs over the Earth history have grown considerably [e.g., 



Bleeker, 



2003 



Harrison et al. 



2005 



Stern, 



2005 



Korenagav 



2006; 



Moizsis et al. 



200 



Van Kranendonk et al. 



2007; 



O 'Neill et al. 



Herz.berg et al. . 



2007 



Condie and Pease . 



2008 



BradleM. 



2008 



Harrison. 



2009; 



2010I1 . Many of previous numerical studies on plate-tectonic convection are, 

27 however, exploratory in nature, and scaling laws relevant to such geological questions are yet 

28 to be established. Given the lack of consensus on why plate tectonics can take place on Eart h 



29 to begin with [e.g., 



2003 



Korenagal 



Moresi and Solomatov . 



1998 



Gurnis et al. . 



2000; 



Bercovici and Karato . 



2007Q, it may be premature to discuss the scaling of plate-tectonic convection, 

31 but it is nonetheless important to seek a strategy to bridge geology and geodynamics by taking 

32 into account peculiar complications associated with plate tectonics. 

33 In this study, I attempt to deri ve the scaling of plate-tecto nic convection using the 



1998], in which the strength of 



34 so-called pseudoplastic rheology nMoresi and Solomatov . 

35 plates is controlled by temperature-dependent viscosity as well as brittle failure. It is known 

36 that, for this approach to be successful, the friction coefficient for brittle deformation has 

37 to be at least one order of magnitude lower than suggested by laboratory experiments. The 

38 presence of pore fluid deep in the oceanic lithosphere is required to explain such low friction, 

39 and because ocea nic lithosphere is like l y to be very dry upon its formation by melting under 



mid-ocean ridges mirth and Kohlstedt 



1996 



Evans et al. 



2005], it may appear to be difficult 



to justify the pseudoplastic approach. The upper half of oceanic lithosphere, however, can 
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42 be pervasively fractured by thermal crackin g, and in the prese nce of surface water, the deep 



200711 . In this mechanism, the strong 



43 hydration of oceanic lithosphere is possible nKorenagd 

44 temperature dependency of mantle rheology actually enhances thermal cracking. Another 

45 concern with the pseudoplastic rheology is that it is determined only by the instantaneous 

46 stress state and does not have any memory to simulate preexisting weakness, though this 

47 limitation is not as grave as it may appear. With the thermal cracking hypothesis, the stiffest 

48 part of oceanic lithosphere is continually damaged as it ages, so preexisting weakness is 

49 globally distributed. Also, whatever the actual weakening mechanism would be, oceanic 
so lithosphere is eventually subducted (on the time scale of 100 Myr), and its memory of 

51 weakness would keep being lost in the deep mantle. For the evolution of oceanic lithosphere, 

52 therefore, the difference between insta ntaneous rheology and history-dependent rheology is 



53 not expected to be vital uackley , 



2000]. In plate tectonics, convective heat loss is dominated 



54 by that from oceanic plates, so even with simple pseudoplastic rheology, we may still hope to 

55 capture the gross characteristics of mantle convection relevant to the long-term evolution of 

56 Earth. 

57 The purpose of this paper is two-fold. First, I will investigate the scaling of plate-tectonic 

58 convection with 'standard' pseudoplastic rheology, which is controlled by friction-based 

59 yield stress and temperature-dependent viscosity. Though there exist a number of numerical 

60 studies using this rheology, the temperature dependency of mantle viscosity is f airly weak in 



61 most of these studies Te.g., 



2004 



Mo resi and SolomatoA 1 1998 



O 'Neill and Lenardid 



Lenardic et all 



2004; 



Stein et al. . 



2007H : the maximum viscosity variation due to temperature 

63 dependency is usually 10 6 . This may not seem to be low because it is high enough to put 

64 convection in the stagnant-lid regime without pseudoplastic rheology. In basally-heated 

65 convection, which is commonly adopted by those previous studies, there is an important 
ee difference regarding thermal structure between stagnant-lid convection and plate-tectonic 

67 convection. In the stagnant-lid regime, most of temperature variations are taken up by the 

68 top thermal boundary layer, so viscosity variation across the top boundary layer is close to 

69 the maximum viscosity variation employed. In the plate-tectonic regime with basal heating, 

70 the top and bottom thermal boundary layers have similar temperature contrasts, i.e., the 
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71 temperature contrast across the top boundary layer (or plates) is basically halved, with the 

72 corresponding viscosity variation of only 10 3 . As explained in more detail later ( §2.11) . the 

73 viscosity contrast across oceanic lithosphere due to temperature dependency is expected to be 

74 at least exp(20) ~5x 10 8 . It is important to use strongly temperature-dependent viscosity so 

75 that we can discuss the scaling of plate-tectonic convection with more confidence. The second 

76 objective of this paper is to discuss the effects of mantle melting by adding depth-dependent 

77 viscosity to the standard pseudoplastic rheology. How mantle melting could modify the 

78 scaling of plate tectoni cs has important implications for the thermal evolution of Earth 



79 nKorenaga . 



2003 



2006], but this issue has not been quantified by fully dynamic calculations, 
so This paper is organized as follows. After describing the details of theoretical formulation 

si (©, I will present numerical results, together with scaling analysis to understand the 

82 systematics of model behavior (©. In the discussion section (@, I will briefly explore how 

83 new scaling laws may be used to inter when plate tectonics initiated on Earth and how it 

84 evolved subsequently. Critiques on previous attempts to derive the scaling of plate tectonics 

85 are also provided. 

se 2. Theoretical Formulation 
s? 2.1. Mantle Rheology 

ss For temperature-dependent viscosity, I employ the following linear-exponential form: 

^ = exp[0(l-r)], (1) 

90 where viscosity is normalized by reference viscosity T}q defined at T* = 1. Temperature is 

91 normalized as 

T — T 

t* = (2) 

93 where T s is the surface temperature (~273 K), and AT is the (arbitrary) temperature scale. The 

94 degree of temperature dependency is controlled by the Frank- Kamenetskii parameter 0, which 



6 



95 can be related to the activation energy E as [e.g. ASolomatov and Moresi l2000|l 



96 



e 



EAT 



R(T S + AT) 



(3) 



Karato and Wu. 



199311 



97 where R is the universal gas constant. For E of ~300 kJ mol -1 [e.g., 

98 and AT of ~1300 K, for example, 6 is ~20. 

99 The coldest part of the lithosphere would be very stiff due to this strongly temperature- 

00 dependent viscosity, but it can also deform by brittle failure. In the continuum limit, this brittle 

01 behavior can be modeled by nonlinear effecti ve viscosity that is adjusted to ensure the stresses 



remain bounded by the yield stress envelope wLoresi and SolomatoA 



1998]. The yield stress 



criterion for brittle deformation may be expressed as 



104 



Ty = CQ + llpQgZ, 



(4) 



05 where cq is the cohesive strength, /i is the friction coefficient, po is reference density, g is 

06 gravitational acceleration, and z is depth. Using the length scale D, which is the depth of a 
fluid layer, and the stress scale T}qk/D 2 , where K is thermal diffusivity, the criterion may be 
nondimensionalized as 
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108 



109 



T +% X Z 



(5) 
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where 



and 



* _ c D z 



WogD 3 



Kilo 



(6) 



(7) 



Using the Rayleigh number defined as 



Ra 



apogATD 3 

Kilo 



(8) 
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us where a is thermal expansivity, the criterion can also be expressed as 

r; = T * + y/?az*, (9) 

us where 

r= 7&r- (10) 

120 In this study, the cohesive strength is assumed to be negligibly small compared to the 

121 depth-dependent component, and Tq is set to T* x 10~ 5 . This is a reasonabl e approxima tion 



Byerlee . 



1978I1 



122 given experimental data on rock friction at low hydrostatic pressure [e.g., 

123 and also allows me to focus on the single parameter y. Note that a nonzero cohesive 

124 strength term appearing for experimental data at high confining pressures could arise from a 

125 pressure-dependent friction coefficient with zero cohesive strength. 

126 The nonlinear effective viscosity for the plastic deformation is calculated as 

r?>* = ^, (ID 

y el 

128 where e^j is the second invariant of the (nondimensionalized) strain rate tensor. The transition 

129 between plastic and ductile deformation is handled by using the harmonic mean of the 

130 temperature-dependent viscosity and the above effective viscosity as 



V*= + • (12) 




132 The effective viscosity for plastic deformation rj* is calculated for any deformation, but when 

133 stresses are smaller than the yield stress (i.e., rj* is large), the harmonic mean above will be 

134 dominated by rjj. 

135 The linear-exponential form of temperature-dependent viscosity [equation (OQ)] predicts 

136 much smaller viscosity variation across the entire lithosphere than the more realistic Arrhenius 

137 form, exp(E /RT), but because the above pseudoplastic rheology effectively eliminates 

138 a drastic viscosity increase in the upper half of the lithosphere, the difference between 
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139 the linear-exponential and Arrhenius forms is actually small \iSolomatoA 12004(1 . It is still 

uo important to use the realistic value of 9 as it controls the strength of the lower half of the 

141 litho sphere. 

H2 When considering the effects of mantle melting, I will add depth-dependent viscosity as 

i« r l*=\ * \ s +— I > (13) 

144 where 

{Art for z* < h* 
(14) 
1 for z* > h* 

146 where h* is the thickness of dehydrated mantle and Arj is a viscosity contrast introduced by 

147 dehydration. 

us 2.2. Governing Equations and Heating Mode 

1 49 The nondimensionalized governing equations for thermal convection of an incompressible 

150 fluid consist of the conservation of mass, 

151 V-u* = 0, (15) 

152 the conservation of momentum, 

153 -VP* + V-[r]*(Vu* + Vu* T )}-RaT*e z = 0, (16) 

154 and the conservation of energy, 

dT* 

155 +u*-vr = V 2 T*+H*. (17) 

at 

156 The unit vector pointing downward is denoted by e z . The spatial coordinates are normalized 

157 by the length scale D, and time is normalized by the diffusion time scale, D 2 / K,. Velocity u* 

158 is thus normalized by k/D. Dynamic pressure P* and heat generation H* are normalized by 
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159 t\qk/D 2 and kAT /(poD 2 ), respectively, where k is thermal conductivity. 

160 In this study, I will focus on thermal convection that is purely internally heated, by 

161 using the insulated bottom boundary condition. There will be no thermal boundary layer at 

162 the bottom, simplifying the scaling analysis of n umerical results. T his heating mode is also 



63 appropriate for the majority of the Earth history nKorenagd 



2008a, section 5.1]. At the same 



164 time, a temperature contrast across the fluid layer is not known a priori, so the maximum 

165 temperature, T^ ax , is not guaranteed to be unity. Some a posteriori rescaling is thus necessary. 

166 The Frank- Kamenetskii parameter is recalculated from its original value 0q as 

e = eoCx 5 (is) 

168 so that exp(0) corresponds to the actual maximum viscosity variation due to temperature 

169 dependency. The internal Rayleigh number may also be defined with 7^ ax as 

Rat = i?ar m * ax exp[0 o (7;; a x " 1)], (19) 

171 in which the total temperature contrast is T£ ax AT, and the internal viscosity, T] ; , is assumed 

172 to be r/oexp[0o(l — ^max)]- Because of purely internal heating, the surface heat flux q is, at a 

173 statistical equilibrium, equal to total heat generation in the fluid divided by surface area 

174 q = p DH, (20) 

175 and the corresponding Nusselt number is calculated as 

q H* 

176 Nu = - — — = • (21) 

kT* ax AT/D T*^ 



The internal heating ratio (IHR), is the difference between heat flux out of the top 



78 bound ary and that into the bottom boundary, normalized by the former [e.g., 



McKenz.ie et all 



1974], i.e., 



Nu top 



10 



182 



183 



and because the bottom boundary is insulated in this study (iVwbot=0), IHR is unity for all runs 
as long as H* is positive. The internal heat production H* does not directly correspond to 
the amount of radiogenic heat production in the mantle, which may be r eferred here as H* ,. 



184 Over the Earth histor y, the mantle has been (usually) cooling with time [Abbott et al. 



Herzberg et al. . 



1994; 



190 



205 



2010], and in the study of mantle convection, this secular cooling is often 
included as part of 'internal' heating. So H* represents both radiogenic heat production and 
secular cooling. Secular cooling is a transient phenomenon, and directly simulating it requires 
us to assume an initial condition for subsolidus mantle convection on Earth, which is hardly 
known. Numerical models for mantle convection are therefore typically run for a number of 
convective overturns to reach a statistical equilibrium so that model results do not strongly 
depend on employed initial conditions. This steady- state modeling approach has to include 
secular cooling as part of internal heating, in order to simulate an Earth-like IHR. The thermal 
evolution of Earth can be studied r easonably w ell by assuming that the mantle is in a quasi 



steady state at each time step re.g.. [DaM 



1980]. 



It is important to distinguish IHR from the convective Urey rati o, U r, which i s the r atio 



of radiogenic heat production in the mantle over the mantle heat flux \\Christensem. 



1985], i.e., 



Ur 



TT* 



Nu 



(23) 



top 



198 The Urey ratio is directly related to the chemical composition of Earth's mantle, and it is a 

199 key parameter to describe the thermal budget of Earth. When Ur is discu ssed, radiogenic heat 



Korenag ca |2008a], 



200 production and secular cooling are considered separately. As noted by 

201 there has been some misunderstanding in the literature by confusing U r with IHR or by 

202 underestimating the sign ificance of secular cool ing, and unfortuately, such confusion still 



seems t o continue Je.g., 



ratio as nKorenagal 



Deschamps et all 



2008a] 



20101 . IHR can be related to the convective Urey 
C c 



Cm ~t~ C c 



l-Ur) 



(24) 



206 where C m and C c are, respectively, the heat capacities of the mantle and the core. The 



11 



209 



207 present-day Urey ratio is probably ~0.2 lKorenagaL\200&d\ . but because the core heat capacity 

208 is only ~l/ 5 of the whole Ear th value, the present-day IHR for Earth's mantle is estimated 



to be ~0.9 nKorenaga 



10 have been higher in the past nKorenaga 



20Q8ah. Based on thermal history considerat ions, the Urey ratio may 



2006 



Herzberg et all 



2010h . so IHR is likely to have 



211 been closer to unity than at present. To first order, therefore, the use of purely internal heating 

212 ((§=1) appears to be a reasonable simplification. 



213 2.3. Notes on Modeling Strategy 

214 Besides the use of pseudoplastic rheology, the numerical model of mantle convection as 

215 specified in the previous sections is kept simple to facilitate the interpretation of modeling 

216 results, and the potential significance of realistic complications, which are neglected in this 

217 study, are discussed in the following. 

218 Because of the insulating boundary condition, bottom heat flux is zero, so there are no 

219 upwelling plumes in the model. The influence of plumes on plate dynamics thus cannot be 

220 examined. Because of the nearly unity IHR expected for Earth's mantle ( §2.21) . however, 

221 such influence may not be of first o rder. The governing e quations employed are based on 

222 the Boussinesq approximation re.g.. \Schubert et all 1200111 . so adiabatic gradients are zero 

223 (i.e., the total temperature contrast 7^* ax Ar is the superadiabatic temperature contrast), and 

224 the model temperature corresponds to potential temperature. The effects of compressibility 

225 on the gross characteristics of thermal convection have been known to be rather minor 



226 Uarvis and McKenzie . 



1980; 



Bercovici et al. . 



199211 



227 For the ductile deformation of the mantle, the Newtonian rheology with linear-exponential 

228 temperature dependency is adopted [equation (Q])], but mantle rheology is known to be much 

229 more complex depending on , at least, stress, pressure, grain size, and chemical composition 



230 [e.g. 



Karato and Wu . 



199311 . In case of pseudoplastic rheology, the difference between the 



231 Arrhenius rheology and its linear-exponential approximation is not important as already 

232 mentioned, and I choose to use the latter because it is specified by only one nondimensional 

233 parameter 6, wh ereas the Arrhenius-type temperature dependency requires three [e.g., 



Korenagal 



3i. 



2009]. Non-Newtonian, stress-dependent rheology can be approximated by 
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235 Newtonian rheology if the activation energy is properly scaled \ Christensen . 



1984]. The 



236 importance of pressure dependence caused by the activation volume is not clear at the 



237 moment. First of all, activation volumes for mant 



238 for upper mantle minerals \lKo renaga and Karato . 



e rhe ology are still poorly known even 



2008]. Second, viscosity increase with 



239 increasing pressure should be at least partly cancelled by viscosity decrease with increasing 

240 temperature along the mantle adiabat. With the Boussinesq approximation employed here, 

241 the use of pressure-independent rheology actually requires non-zero activation volume, the 

242 effect of which is assumed to be cancelled exactly by temperature v ariations alo ng the adiabat. 



Solomatov 



1996], but 



243 Grain size variation can affect mantle dynamics considerably [e.g., 

244 how grain size should evolve in the convecting mantle is still poorly understood, so it appears 

245 premature to consider its effect in this study. The effect of composition on mantle rheology 

246 is taken into account when dehydration stiffening is effected by dept h-dependent vi scosity 



247 [equation (1141)1 . There are of course other compositional effects [e.g., 



248 effect of dehydration appears to be most important at l east for the upper mant 

,0; 



[e.g., 



Karato et al. 



1986; 



Mei and Kohlstedt 



2000a, 



Karatd 1200811 . but the 
e rheology 



Faul and Jackson 



2007], and mantle 



1996]. 



250 dehydration is always expected whenever mantle melts jHirth and KohlstedA 

251 Another important rheological aspect for large-scale mantle dynamics is a viscosity jump 

252 at the base of the upper mantle, which has been estimated to be on t he order of ^ 10-100 



253 primarily through the geodynamical modeling of Earth's geoid [e.g., 



Hager , 



198411 . Su ch 



254 inference is, however, a lso known to suffer from considerable nonuniqueness [e.g., 



Kido and Ca dek. 



King . 



1995 



199711 . and the viscosities of the upper an d lower mant l e may not be very 



2009]. Furthermore, 



256 different if the mantle transition zone has a lower viscosity nSoldati et all 

257 even if the lower mantle does have a higher viscosity than the upper mantle, it applies only 

258 for the present-day situation. When the mantle was hotter in the past, the viscosity contrast 

259 may be smaller or even reversed if the lower mantle rheology is more temperature-dependent 

260 (i.e., higher activation energy) than the upper mantle counterpart. Rheological stratification in 

261 the mantle is an important subject, but these uncertainties imply a variety of situations to be 

262 considered, so it is left for future studies. 



The mantle transition zone is also characterized by multiple phase transitions, and 
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264 in particular, the effects of the endothermic phase change at the base of the transition 



266 



267 



Chri stensen and Yuen , 



1984; 



Tackley et all 



1993; 


Solheim and Peltier , 


1994; 


Yuen et al. 



1994]. Numerical studies with strong plates exhibit, however, only a mo dest influence of 



268 endothermic phase change on mantle dynamics [e.g., 



Zhong and Gurnis 



1994J, and recent 



269 experimental studies further suggest that the Clapeyron s l ope of the endot hermic phase change 



is likely to be only —1.3 MPa K nKatsura et al. 



2003; 



Fei et al. . 



2004], which is much less 



271 negative than previously thought. Modeling phase transitions, therefore, is not considered to 

272 be essential. 

273 Finally, the model is 2-D Cartesian, whereas the use of a 3-D spherical shell would be 

274 most appropriate. The restriction to 2-D modeling is primarily to generate a large number 

275 of modeling results (with modest computational resources) so that scaling analysis becomes 

276 more robust, though I do not expect scaling laws to change dr astically by mov i ng fro m 2-D to 



277 3-D. The effect of sphericity would likely be of minor nature nBercovici et al. 



278 on iso viscous convection models using 3-D spherical shells, for example, \Deschamps et al. 



200 0[| . Based 



279 Il2010ll recently derived the following heat-flow scaling (adapted here for the case of purely 

280 internal heating), 

Nu « 0.59 / o.05 jRa o.300-o.oo 3/) (25) 

282 where / is the ratio of the core radius to the total radius of a planet. The ratio / is 0.55 for 

283 Earth and unity for Cartesian, so it can be seen that sphericity has virtually no impact on this 

284 scaling. 

285 The convection model of this study is, therefore, simple but probably not simpler than 

286 necessary. In any event, this study should provide a reference point, by which the effects of 

287 any additional complication can be quantified in future. 



288 

289 
290 



3. Numerical Results and Scaling Analysis 

The finite element code of \Korenaga and Jordan [1200311 was used to solve the coupled 
Stokes flow and thermal advection-diffusion equations (fT5T)-(fT7l). The benchmark tests of this 
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291 code c an be found 'm \Korenaga and Jordan Q2003D for Newtonian rheology and in \Korenaga 

292 j2009l for non-Newtonian rheology. To reduce wall effects, the aspect ratio of the convection 

293 model is set to 8, and the model domain is discretized with 400 x 50 uniform 2-D quadrilateral 

294 elements. With this mesh resolution, model parameters are chosen so that Nu does not 

295 exceed 20 and the top thermal boundary layer contains at least a few elements vertically on 

296 average. The nondimensional surface temperature is fixed to zero, and the bottom boundary is 

297 insulated. The top and bottom boundaries are free-slip, and a reflecting boundary condition is 

298 applied to the side boundaries. In all cases, Ra is set to 10 , but Rai varies greatly because of 

299 different combinations of 0q and H* (and thus T^L). 

300 The initial temperature condition is specified as 



301 T*(x\z*) =z* + acos(7W*)sin(7z;z*) +e, (26) 

302 where a is usually 0.2, and e is random fluctuation with the amplitude of 10~ 3 . When the 

303 assumed mantle rheology is appropriate for the operation of plate-tectonic convection, this 

304 initial condition quickly brings the system to that mode of convection. Otherwise, the system 

305 gradually migrates into the mode of stagnant-lid convection. If I start with a uniformly hot 

306 fluid instead, the system always begins with stagnant-lid convection, and sublithospheric 

307 mantle is heated up considerably until the onset of plate tectonics. Very low viscosity (and 

308 thus very high convective velocity) beneath the stagnant lid during this initial period means 

309 exceedingly small time steps for numerical integration, so this type of initial condition is 

310 not computationally efficient when aiming at statistically steady states required for scaling 

311 analysis. 

312 3.1. Convection Diagnostics 

313 A typical snapshot of model run is shown in Figure [T^. This is the case of y=0.6, 



314 0o=15, and H*=20, with the standard pseudoplastic rheology [equation (fl2)) 1. In addition to 

315 Zj^x, I calculate two more measures for the temperature scale. One is the domain-average 
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316 temperature, 



317 



320 



324 



[T*} = jj T*dx*dz* I jjdx*dz* 



(27) 



318 and the other is the (self-consistent) internal temperature nKorenagaL I2009M . 



T* = YZJijg, I / T * dx * I I dx * I dz * 



(28) 



where 8' = T*/H*. 

321 The vigor of convection can be quantified by calculating the root-mean- square 

322 velocity vj^, and a velocity diagnostic most indicative of the mode of convection is the 

323 root-mean-square surface v elocity v*. To quantify how plate-like the surface velocity field is, 



Weinstein and Olson n 199211 introduced the notion of 'plateness', and for convection exhibiting 



325 multiple plates with different velocities, I use the following definition of plateness, 



dx* 



e'<x 



dx* 



(29) 



327 where 



328 



337 



dv*(z* = 0) 



dx* 



(30) 



329 The parameter P x measures the fraction of surface with normalized strain rate e' smaller than 

330 the given threshold x. The velocity profile shown in Figure [TJ5, for example, has Pq.i of 0.76 

331 (Figure [Tt). Like other definitions of plateness, P x varies from to 1, with higher values 

332 corresponding to more rigid behavior. For comparison, actual plates on Earth tend to have 

333 wide diffuse boundary zo nes, which occupy ~15% of the surface area at the present day 



334 \Gordon and Stein . 



1992]. 



335 To understand the spatial distribution of viscous dissipation, I define 4> c . as the viscous 

336 dissipation within the region above z* = c, 



[ Q^e^dx^dz* 



(31) 



338 and calculate 4>^, $5 n, and <I> (= 4>i), where 5 is (on average) the maximum thickness of the 
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339 top thermal boundary layer [e.g., [S^Ma LL967|], 



340 



8 = 2Nu~\ (32) 



341 Here 5 is nondimensionalized by the model depth D. 

342 Each case was run up to t* = 6. The cumulative heat generation and cumulative heat loss 

343 from the surface are monitored, and when these two start to match within ~1%, I judge that 

344 the system has reached statistically steady state. This usually takes place at t* =~ 2 — 3, and I 

345 use subsequent model results to calculate time-average values of key diagnostics such as Nu 

346 and 4>. The (one) standard deviation of time-averaged Nu is typically less than 1%, whereas 

347 the standard deviation of Nu itself is often greater by one order of magnitude, reflecting the 

348 highly time-dependent nature of convection (Figured]!). Surface velocity exhibits even greater 

349 time dependency (Figure QJ). 

350 3.2. Reference Scaling 

351 A total of 82 cases were run with the standard pseudoplastic rheology using different 

352 combinations of 7 (0.1-1), Qq (10-25), and H* (8-20). The summary of convection diagnostics 

353 is reported in Tables IA1I and IA21 eight runs resulted in stagnant-lid convection, and others 

354 exhibited plate-tectonic convection. 

355 Three different temperature scales, 7^ ax , (T*), and T*, are correlated well to each other 

356 (Figured). Regardless of the mode of convection, 7^ is distinctly higher than (T*), and this 

357 is because the thickness of the top thermal boundary layer is not trivial in those runs. The 

358 maximum Nu achieved is only ~20 (Figure Wp), and these different temperature scales are 

359 expected to converge as Nu increases. Stagnant-lid runs are characterized by similar 7^ ax and 

360 T*, because of small temperature variations beneath the stagnant lid. 

361 The relation between Nu and Rat appears to roughly follow the classical scaling of 

1 /3 

362 Nu °z Ra i 1 within runs with the same 7 and similar (Figure [2b)- Varying 7 has considerable 

363 effects on the scaling of Nu as well as v* ms (Figure [2fc). The distinction between plate-tectonic 

364 and stagnant-lid runs is very clear in the correlation (or lack thereof) between and v* 
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365 (Figure O. Average surface velocity in these plate-tectonic runs is higher than corresponding 

366 v* ms because the latter involves averaging over the entire domain, the majority of which moves 

367 more slowly than surface plates. 

368 My scaling analysis to understand the syste matics of these model results is based on 



373 



1966TI . Because of pseudoplastic 



369 the local stability of top thermal boundary layer ftHowardl 

370 rheology, the effective viscosity of the top boundary layer or the effective lithospheric 

371 viscosity, r\i, is expected to be higher than the interior viscosity, and I denote the viscosity 

372 contrast between them as 



At]l = TlL/rii. (33) 



374 Viscosity in the top thermal boundary layer varies considerably as specified by equation (PT2l) . 

375 and the effective lithospheric viscosity is an attempt to capture the overall stiffness of the 

376 boundary layer by just one viscosity value. For the stiff boundary layer to subduct, it has to 

377 become convectively unstable at least, and by assuming that the maximum thickness of the 

378 boundary layer 8 D corresponds to marginal stability, the following relation should hold: 

q W (C^-)(gP) 3 

380 where Ra c is the critical Rayleigh number. By using the relation between 8 and Nu 

381 [equation (|32l) l. this marginal stability criterion may be rearranged as 

D \ 1/3 

RaA ' . -i/3 



382 Nu = 2(— -j AT7 L 1/J , (35) 



383 or 

A7lL = R^- (36) 

385 Hereinafter Ra c is set to 10 3 . 

386 Equation d36l) may be regarded as a way to extract lithospheric viscosity contrasts 

387 from the measured pairs of Rat and Nu. The lithospheric viscosity contrast calculated this 

388 way increases as 9 increases, and this 6 sensitivity is greater for higher 7 (Figure [3^). The 
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389 following functionality appears to be sufficient to reproduce the first-order behavior of the 

390 viscosity contrast, 

ATj L (y,e)=exp[a(y)0], (37) 

392 which converges to unity at the limit of zero 6. The coefficient a(y) is determined by linear 

393 regression for each group of runs with the same y (Figured). Excluding the result for y of 

394 0.1, the runs with which are characterized by rather low plateness (Table lAll) . the coefficient 

395 is linearly correlated with y in the logarithmic space (Figure[3t>), which may be expressed as 



a{y) ^0.327y°- 647 . (38) 



397 Equations (|37l) and 081) can predict At]^ reasonably well over the range of four orders of 

398 magnitude (Figure [3]:). The prediction for Nu through equation (1351) has the average error 

399 of ~10% (Figure [3ti). This error is considerably larger than that observed for the heat-flow 



400 scaling derived by 



Moresi and Solomatow II 199811 . which was based on convection with the 



401 aspect ratio of one. The use of the wide aspect ratio (= 8) in this study and resultant time 

402 dependency in convection patterns may be the source of these scatters. 

403 The lithospheric viscosity contrast At/l calculated from equation (l36l) exhibits broad 

404 correlations with other convection diagnostics (Figure H]). Higher Arji generally gives rise 

405 to higher plateness (Figure Hk) and greater viscous dissipation in the top boundary layer 

406 (Figure Ht). How viscous dissipation is distributed within the boundary layer, however, seems 

407 to be insensitive to variations in At/l as ~ 0.5-0.6 for all of plate-tectonic runs 

408 (Figure |4]1). 

409 Assuming the half-space cooling of lithosphere, the (average) maximum plate thickness 

410 8 is related to the average length of plates, L, as 

1/2 



s 



412 where Lj (Dv*) is the average time from a ridge to a subduction zone. Thus, the average aspect 
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13 ratio of convection cells may be calculated from Nu and v* as 



418 



D Nu 2 ' 



(40) 



415 for which equation (|32l) is used. The aspect ratio gradually increases as Arji increases, i.e., 

416 stronger plates tend to be longer (Figure @})). By assuming some empirical relation for L/D 

417 (e.g., L/D ~ AriV 6 ), equation (l40l) may be rearranged as scaling for v*, 



D \Ra, 



^\ 2/3 Ar?L 2/3_ 



(41) 



419 That is, unlike the scaling for heat flux [equation (|35l) l. some information on the aspect ratio 

420 of convection cells is essential for the scaling for surface velocity. 



421 3.3. Effect of Shallow Stiffening 

422 On Earth, the creation of new plates at mid-ocean ridges is usually accompanied by the 

423 melting of upwelling mantle (unless the mantle is too cold), and this chemical differentiation 

424 along the gl obal mid-oce an ridge system constitutes the dominant fraction of terrestrial 



425 magmatism nCrispl 



1984]. This mantle melting results in the formation of oceanic crust as 



426 well as depleted mantle lithosph ere, both of which are chemical 



427 to the underlying asthenosphere nOxburgh and Parmentier 



y more buoyant with respect 



197711 . and the depleted lithosphere 



428 also becomes intrinsically m ore viscous (by ~10 3 ) because of dehydration caused by melting 



429 \Hirth and Kohlstedt 



1996|]. As long as plate tectonics is taking place, the chemical buoyancy 



430 of oceanic lithosphere is insignificant as resistance to subduction because the basalt-to-eclogite 

431 transition at relatively shallow depth (< 60 km) makes the subducting slab compositionally 



432 denser than the surrounding mantle \e.g.. \Ringwood and Irifune 



19880. In this study, therefore, 



433 I fo cus on the effec t of dehydration stiffening on the scaling of plate-tectonic convection. As 



434 in 



Korenaga 



[2009], instead of tracing the advection of the dehydrated slab through time, I 

435 use the depth-dependent viscosity that is fixed in time [equation (fT4l)1. Though the subducting 

436 slab loses the extra viscosity contrast At] as soon as it passes the given depth h*, the effect of 
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437 shallow stiffening on slab bending can still be evaluated with this scheme. 

438 The scaling of Nu [equation (1351) 1 indicates that different combinations of Rai and Ar\i 



439 can produce the same Nu. Thus, in terms of the efficiency of heat transport, a run with high 

440 Rat and high At]^ may be indistinguishable from that with low Rat and low Arj^, but shallow 

441 stiffening may affect these cases differently. A wide variety of plate-tectonic cases were 

442 thus simulated by varying y (0.4-0.8), 0q (10-25), and H* (2-20), and shallow stiffening was 

443 incorporated with h* ranging from 0.1 to 0.3 and Arj ranging from 3 to 10 3 . Stiffening by 

444 mantle melting is limited mostly to the top 200 km or so (i.e., h* < 0.07), and the use of 

445 greater h* is to study the asymptotic behavior of stiffening effects. Also, the top thermal 

446 boundary layer in numerical modeling is thicker than actual oceanic lithosphere because of 

447 relatively low Rat used in this study, so h* has to be comparably large in order to reproduce an 

448 Earth-like combination of thermal and compositional boundary layers. The number of runs is 

449 225 in total, with 19 stagnant-lid runs (Tables IA31A8I) . Convection diagnostics were measured 

450 in the same way for reference runs. 

451 The scaling analysis for runs with shallow stiffening is founded on that for the reference 

452 runs. As in the previous section, the lithospheric viscosity contrast ArjL is calculated from the 

453 measured pair of Nu and Raj. A key issue is how this viscosity contrast is influenced by the 

454 additional depth-dependent viscosity, and this influence may be measured by the deviation 

455 from the prediction based on equation (1371) . The predicted viscosity contrast is based solely on 

456 y and 6, and it is denoted as A7|L, ref to distinguish from the actual Arji. The ratio Arji/ Arj^ TC f 

457 is loosely correlated with Ar] and h* as one may expect (Figure [5^,b); higher Ar\ or h* leads to 

458 higher Ar\i than predicted by equation (1371) . A better correlation may be seen between the two 

459 ratios, Ar\i/ Ar\i mf and h*/h* ef (Figure[5t), where h* ef is defined as 



460 



h 



ref 



.* 



Nu, 



-l 

ref J 



(42) 



461 and 



462 




(43) 
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463 The parameter h* ei is the averaged thickness of the top thermal boundary layer expected for 

464 a run with the same Rcii but with the standard pseudoplastic rheology. The ratio Arj^/ 'Arj£ ref 

465 increases as h* //z* ef but eventually saturates and never exceeds the given At]. This behavior 

466 may be represented by the following functionality, 



468 which means that Ar\^ converges to the simple product of Ar]^ ref and Ati when h* is 

469 sufficiently greater than h* ef . The parameter x controls how high h* should be with respect 

470 to h* ef in order to achieve the convergence; greater x means that thicker h* is required. By 

471 trying a range of values, I found that x ~ 6 can reproduce the measured Ar\ L reasonably well 

472 (Figure (5]i). 

473 The relations between Arji and other convection diagnostics are more ambiguous than 

474 observed for the reference runs (Figure©. It appears to be premature to parameterize the 

475 aspect ratio L/D as a simple function of the lithospheric viscosity contrast (Figure [6})), and 

476 more thorough work is clearly required to better understand the self-organization of plate 

477 tectonics. It is still interesting to note, however, that high plateness is possible even with low 

478 viscosity contrast (Figured) and that plate-tectonic convection can occur even when most of 

479 viscous dissipation takes place in the top boundary layer (Figure©). 

480 3.4. Conditions for Plate-Tectonic Convection 

481 The condition for plate-tectonic convection is found to be seen most clearly in the 

482 covariation of Rat and Ar\i (Figure [7]). Plate-tectonic convection is possible even with high 

483 lithospheric viscosity contrast if Raj is sufficiently high, and the critical viscosity contrast, 

484 above which plate-tectonic convection is unlikely, appears to be 



467 




(44) 



Al]L )C rit ~ 0.25i?a. 



1/2 



(45) 



485 



486 though there are a few runs that slightly violate this threshold. 
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487 To understand the meaning of this scaling, I consider stress balance at the bending of 

488 subducting slab, which is probably the most critical part of plate-tectonic convection. First, 

489 the stress due to the negative buoyancy of the slab may be expressed as 

T S ~ apog(T* ax AT)D, (46) 

491 or by normalizing the (internal) stress scale, rjjK/D 2 , 

492 t| ~ Ra t . (47) 

493 Second, the bendi ng stress should be propor tional to the lithospheric viscosity and bending 



494 strain rate as [e.g., 



Conrad and Hagen 



1999I1 



v s (8D) 

t b ~ Il^t 1 , (48) 
496 where R is the radius of curvature, and its nondimensionalized form is 

tf~Aifcv;5QQ ocA^X 1/3 (^) • (49) 

498 Finally, by assuming t s rs t b at At]l = A7]L, crit , we may derive the following scaling for the 

499 radius of curvature, 

500 -ocR a 7 l/6 . (50) 

D 

501 Thus, the radius is weakly dependent of the vigor of convection, and it becomes smaller 

502 for more vigorous convection. Without this variation in the radius of curvature, the critical 

503 viscosity contrast would be more sensitive to a change in Rat (i.e., proportional to Rat instead 

504 of Ra 1 / 2 ). 

505 Note that all of 307 runs reported here are either strictly plate-tectonic or stagnant-lid 

506 convection, and there is no case of episodic overturn, in which the system periodically 



507 goes back and forth between plate-tectonic and stagnant-lid modes wLoresi and SolomatoA 
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511 



513 



525 



508 LL998]. This is consistent with the use of virtually zero cohesion strength and finite friction 

509 coefficient in this study ( §2.11) . because the possibility of the episodic overturn mode appears 



510 to be important only with nontrivial c ohesion strength [e.g., 



Stein et al. . 


2004; 


O'Neill et al. 


2007 1 



Moresi and Solomatov 



1998 



2007]. The use of purely internal heating (thus the lack of 



upwelling plumes) in this study might also be responsible. 

4. Discussion and Conclusion 



514 On the basis of the scaling of Nu [equation ([351) 1 and the parameterization of At]l 

515 [equations (|37l) . (1381) . and (l44l)l. it is now possible to discuss the relation between mantle 
sis temperature and surface heat flux, which is fundamental to our theoretical understanding of 
517 the long-ter m evolution of Earth. Becaus e some of key model parameters are still poorly 



sis known [e.g., 



Korenaga and Karato 



2008], the following exercise should be regarded as a 

519 preliminary case study. As explained below, the self-consistent construction of a plausible 

520 heat-flow scaling law requires modeling the thermal and chemical evolution of Earth at the 

521 same time, so a more extensive exploration of the scaling of plate tectonics will be reported 

522 elsewhere. 

523 First, for the dependency of viscosity on mantle potential temperature T p , the following 



524 Arrhenius form is used, 



Vt(T p ) 



( E 
T] r exp I 



RT n 



E 

RT r 



(51) 



526 where r\ r is reference viscosity at T„ = T r , and E is assumed to be 300 kJ mol nKorenagal 



2006]. The reference temperature is set to 1623 K (1350' 



C), which corr e spond s to the 



2007]. The 



528 present-day potential temperature of the ambient mantle Werzberg et all 

529 Frank-Kamenetskii parameter 6 is calculated from equation ©, and the reference lithospheric 

530 viscosity contrast Ar^^f is calculated with 7 of 0.8, which corresponds to the effective friction 

531 coefficient of ~0.02 [equation (fT0l)1. 

532 The internal Rayleigh number Raj is then calculated with the above temperature- 

533 viscosity and the following values: a=2 x 10~ 5 K _1 , po=4000 kg m~ 3 , g=9.8 m s -2 , and 

534 D=2900 x 10 3 m. The dehydration of the mantle beneath mid-ocean ridges is assumed to 
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535 take place when the upwelling mantle crosses the solidus for dry pyrolitic mantle, and the 

536 i nitial pressure of melt ing (in GPa) can be calculated from the potential temperature (in K) as 



538 



544 



537 IKorenaga et all 



2002] 



1423)/100, 



(52) 



539 and the thickness of dehydrated mantle h m is given by P /(p g). The nondimensional 

540 thickness h* is h m /D, and the viscosity contrast due to dehydration A77 is assumed to be 10 2 

541 here. The lithospheric viscosity contrast At]l gradually increases for higher T p (Figure [8fc, 

542 case 1) because hotter mantle starts to melt deeper (Figure [8b), but the effect of shallow 

543 stiffening on the viscosity contrast saturates at T p ~ 1600°C, above which the contrast slightly 
decreases because of smaller 6 for higher T p [equation ©]. 

Lastly, the global heat flux Q is calculated as 



546 



560 



Q = JcATpNu/D, 



(53) 



547 where A is the surface area of Earth, and k is assumed to be 4 W m 1 K . The reference 

548 viscosity r\ r is set to 10 19 Pa s so that the predicted global heat flux matches the present-day 



549 convective heat flux of ^38 TW jKorenag a. 



2008ah (Figure [8]1, case 1). The effect of shallow 



550 stiffening suppresses the heat flux considerably and even reverts the sense of temperature 

551 sensitivity; the flux is lower for higher T p above ~1450°C until the effect of shallow stiffening 

552 becomes saturated at ~1600°C. For the temperature range of 1 350-1600°C, which is most 



553 relevant to the thermal evolution of Earth for the last 3.5 Gyr jHerzberg et al. . 



2010], the 



554 predicted rel ation between the mantle temperature and surface heat flux closely resembles that 



555 suggested by 



Korenag a 



Q2OO60 along a similar line of reasoning but on the basis of the global 

556 energy balance. 

557 Note that the reference viscosity r\ r of 10 19 Pa s (at present-day potential temperature) 

558 may be appropriate for asthenosphere but would typically be regarded as too low to represent 

559 the whole mantle, for which the viscosity of 10 21 -10 22 Pa s is usually assumed. The effective 



lithospheric viscosity contrast is ~3x 10 2 in this example (Figure^), which is comparable to 



25 



561 the discrepancy. Traditionally, the average viscosity of the present-day mantle is estimated to 

562 be of that magnitude, in order to explain surface heat flux (or equivalently, plate velocities) 



ses [e.g., 



Hager 



1991 



Bercovici et al. 



2000], because the aforementioned geoid-based studies 



564 ( §2.31) can constrain only relative variations in viscosity and are insensitive to the absolute 

565 values of viscosity. The heat-flow scaling of equation (|35l) suggests that a lithospheric 

see viscosity contrast alone could regulate surface heat flux without invoking a viscosity increase 

567 in the lower mantle. 

ses The use of constant viscosity contrast for dehydration stiffening At~i for the entire 

569 temperature range (thus implicitly over the entire Earth history) is equivalent to assuming 

570 that the water content of the convecting mantle does not change with time. If the mantle is 

571 drier than present, for example, the viscosity contrast would be smaller, and if the mantle is 

572 completely dry, mantle melting should not cause any viscosity change. By combining the 

573 thermal budget of E arth with geolog ical constraints on sea level change and with the growth 



574 of continental crust, 



Korenagd [|2008b|l suggested that the volume of Earth's oceans is unlikely 



575 to have been constant with time and that the mantle may have been gradually hydrated by 

576 subduction starting with a very dry state in the Archean. As the second example (denoted 

577 as case 2 in Figure [8]), I consider effective heat-flow scaling expected for this scenario. For 

578 simplicity, the mantle is assumed to have been hydrated linearly from the co mpletely dry state 



579 to the present state, as it cooled from 1550°C, to 1350°C for the last ~3 Gyr nHerzberg et all 

580 |20 1QTI . For T p greater than 1550°C, therefore the internal viscosity is intrinsically higher by 

581 At], and this viscosity contrast gradually diminishes as T p approaches 1350°C. This is reflected 

582 in how Rai varies with T p (Figured, case 2). At the same time, the viscosity contrast due to 

583 mantle melting is unity at T p > 1550°C and gradually increases to the full value Ar] at T p of 

584 1350°C. The total lithospheric viscosity contrast At]l in this scenario is much reduced than the 

585 previous example (Figure[8fc), but because of the overall reduction in Rai, the surface heat flux 
see is suppressed further, and the inverse relationship between mantle temperature and heat flux 

587 dominates heat-flow scaling during the mantle hydration period (Figure (8]i). Obviously, this 

588 type of calculation should be done more self-consistently by modeling the thermal evolution of 

589 Earth together with its global water cycle, and what is presented here is only a crude estimate. 
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590 In both cases, the lithospheric viscosity contrast is always smaller than its threshold 

591 (Figure [St), so plate-tectonic convection seems to be dynamically plausible throughout 

592 the Earth history, as long as surface water ex ists to hydrate the litho sphere and reduce the 



593 



effective friction coefficient \Korenaga . 



2007]. In particular, the gradually hydrating mantle 

594 (case 2) helps to maintain relatively small Arj^ even with deeper mantle melting at higher T p , 

595 facilitating the operation of plate tectonics in the early Earth. 



Though previous attempts to estimate the heat-flow scaling of plate tectonics nKorenagcL 



2003, 



)6] 



2006] have already predicted the inverse relation between mantle temperature and 



598 surface heat flux as indicated by Figure[8l there are a few important differences. First, because 

599 the effect of shallow stiffening eventually saturates [equation ([441) 1. the inverse relation is 

600 restricted to a certain temperature range. This subtle behavior is difficult to derive from the 

601 global energy balance approach adopted by the previous studies. Second, the global energy 

602 balance can be exploited to derive heat-flow scaling by assuming the mode of convection, so 

603 whether plate-tectonic convection is plausible or not cannot be addressed. Finally, the energy 

604 balance approach has a few poorly constrained parameters, such as the radius of curvature 

605 for plate bending, effective lithospheric viscosity, and the aspect ratio of conve ction, and it is 



606 possible to obtain wildly different results by varying them independently \e.g.. \Davies . 



2009]. 



607 Though similarly suffering from parameter uncertainty (e.g., y and Arj) and from the very 

60s assumption of the pseudoplastic rheology as well, the present study provides a fully dynamic 

609 framework in which heat flow, velocity, lithospheric viscosity, aspect ratio, and the radius of 

610 curvature are all connected in a self-consistent manner. 

en Appendix A: Numerical results 

612 Tabulated are selected convection diagnostics for statistically steady-state solutions as 

613 described in the main text. The Rayleigh number Ra is 10 6 for all cases. Reference runs refer 

614 to calculations with the standard pseudoplastic rheology, and runs with shallow stiffening 

615 refers to those with additional depth-dependent viscosity. 
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Figure 1. Example of simulation results from the case of 7=0.6, 0o=15, and H*=20 with the 
standard pseudoplastic rheology [equation (fT2l)l. (a) Snapshot of the temperature field. Darker 
shading corresponds to lower temperature, (b) Surface velocity profile and (c) horizontal strain 
rate scaled by by the average surface velocity, corresponding to the snapshot shown in (a). Pq i 
is 0.76 for this particular velocity profile, i.e., 76% of the surface has the scaled strain rate 
lower than 0.1 (shown by dotted in (c)). (d) Nusselt number and (e) root-mean- square surface 
velocity as a function of time (shown in gray). For this model run, statistically steady state 
was achieved at t* = 1.8, and running average is taken from the subsequent model results. In 
(d) and (e), the running average and its uncertainty (1 o) is shown as solid and dotted lines, 
respectively, and the one standard deviation of the temporal variation itself is shown by dashed 
line. 
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Figure 2. Correlations among convection diagnostics from reference runs, (a) Internal tem- 
perature T* (cross) and maximum temperature 7^ ax (solid circle) are compared with domain- 

1 /3 

average temperature (T*). (b) Nu and Rat. Dashed line indicates the slope of Ra/ . (c) v* ms 
and Rcii. (d) v* and v^. In (b)-(d), different symbols denote runs with different y, as shown by 
the legend in (b). 
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Figure 3. The systematics of reference runs can be summarized through effective lithospheric 
viscosity contrast ArjL. (a) Arji as a function of 6. Different symbols denote runs with different 
7, and dashed lines are fitted trend in the form of equation (1371) for each y. Data with ArjL 
greater than 10 4 are stagnant-lid runs, which are excluded from linear regression, (b) The 
fitted coefficient a as a function of y Dashed line represents equation (1381) . (c) Comparison 
of measured ArjL with predicted values based on equations (|37|) and (1381) . (d) Comparison of 
measured Nu with predicted values based on equation (1331) . 
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Figure 4. Covariations of measured litho spheric viscosity contrast Arji with (a) plateness Pq \, 
(b) the average aspect ratio of convection cells L/D, (c) the fraction of viscous dissipation 
taking place in the top thermal boundary layer with respect to that in the entire domain <I>g /4>, 
and (d) the fraction of viscous dissipation in the upper half of the boundary layer with respect 
to that in the entire boundary layer <&§ n/^s- Dashed line in (b) corresponds to L/D = Arj^ 6 . 
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Figure 5. The effect of additional depth-dependent viscosity on the lithospheric viscosity con- 
trast. The deviation from the prediction based on standard pseudoplastic rheology, At]l/ AT]L ire f, 
is shown as a function of (a) the viscosity contrast due to dehydration At], (b) the thickness of 
dehydrated layer h*, and (c) the same thickness but scaled by the reference thickness, h*/h* ef . 
Different symbols denote different h* in (a) and different At] in (b) and (c). Note that in (a) and 
(b) the values of At] and h* are slightly perturbed randomly for display purposes, (d) Compar- 
ison of measured At]l with predicted values based on equation (|44"1) . for three different values 
of*. 
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Figure 6. Same as Figure|4]but for runs with shallow stiffening. Different symbols correspond 
to different groups of runs: 7 of 0.4 with high H* (solid circle; Tables IA3HA5I) . 7 of 0.4 with 
low H* (gray circle; Tables IA6HA7]), and 7of 0.8 (open circle; Table|A8])- 
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Figure 7. Covariation of i?a, and Arji for all model runs. Solid and open symbols denote 
plate-tectonic and stagnant-lid runs, respectively. Dashed line represents an approximate divide 

1/2 

between these two modes of convection (At]l ~ 0.25i?a- ). 
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Figure 8. A worked example of how the new scaling laws of plate-tectonic convection may 
be used to build heat-flow scaling for Earth, (a) Internal Rayleigh number Rcii as a function of 
mantle potential temperature T p . (b) Thicknesses of dehydrated lithosphere (h m , dotted) and 
reference thermal boundary layer (/z re f, solid for case 1 and dashed for case 2) [equation (|42l) l. 
(c) Lithospheric viscosity contrast for case 1 (thick solid), case 2 (thick dashed), and a reference 
case with no effect of mantle melting, i.e., h* m = (thin solid). Also shown are the critical 
viscosity contrast for plate-tectonic convection (dashed for case 1 and dotted for case 2). (d) 
Relation between T p and surface heat flux Q. Legend is the same as in (c). Star denotes 
convective heat flux at the present day (38 TW at 1350°C). See the main text for details. 
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799 Tables 



Table Al. Numerical results for reference runs (y < 0.5) 



r 


6o 


H* 




7.* 

i 


<n 


Nu 





Rat 




v rms 


v s * 




P0.1 


4> 


<S> S /(S> 


0.1 


15 


8 


0.86 


0.59 


0.60 


9.27 


12.9 


l.lOx 


10 5 


6.40 x 10 1 


1.49 x 


10 2 


0.33 


2.73 x 10 7 


0.19 


0.1 


15 


12 


0.93 


0.68 


0.69 


12.93 


13.9 


3.18x 


10 s 


1.12 x 10 2 


2.26 x 


10 2 


0.30 


4.25 x 10 7 


0.14 


0.1 


15 


16 


0.98 


0.75 


0.76 


16.39 


14.6 


6.85 x 


10 5 


1.68 x 10 2 


3.10 x 


10 2 


0.26 


5.77 x 10 7 


0.13 


0.1 


15 


20 


1.02 


0.82 


0.81 


19.67 


15.3 


1.32 x 


10" 


2.33 x 10 2 


3.97 x 


1 0- 


0.22 


7.31 x 10 7 


0.14 


0.1 


20 


12 


0.96 


0.71 


0.70 


12.50 


19.2 


4.36 x 


I0 5 


1.00 x 10 2 


2.22 x 


I0 2 


0.27 


4.22 x 10 7 


0.16 


0.1 


20 


16 


1.00 


0.77 


0.77 


15.95 


20.1 


1.08 x 


10 6 


1.57 x 10 2 


3.18x 


1 0^ 


0.28 


5.77 x 10 7 


0.13 


0.1 


20 


20 


1.04 


0.84 


0.82 


19.28 


20.8 


2.21 x 


10" 


2.19 x 10 2 


4.02 x 


in 2 


0.25 


7.30 x 10 7 


0.13 


0.1 


25 


12 


0.97 


0.74 


0.70 


12.32 


24.4 


5.15 x 


10 5 


8.78 x 10 1 


2.13 x 


10 2 


0.22 


4.19 x 10 7 


0.16 


0.1 


25 


16 


1.02 


0.80 


0.77 


15.69 


25.5 


1.68 x 


10 6 


1.43 x 10 2 


3.01 x 


10 2 


0.27 


5.73 x 10 7 


0.14 


0.1 


25 


20 


1.05 


0.85 


0.82 


19.04 


26.3 


3.74 x 


10 6 


2.06 x 10 2 


3.97 x 


10 2 


0.28 


7.28 x 10 7 


0.13 


0.2 


15 


8 


0.94 


0.66 


0.66 


8.48 


14.2 


4.05 x 


10 5 


5.97 x 10 1 


1.11 X 


10 2 


0.45 


2.67 x 10 7 


0.23 


0.2 


15 


12 


1.00 


0.74 


0.75 


12.04 


15.0 


9.57 x 


10 5 


1.21 x 10 2 


2.06 x 


10 2 


0.48 


4.20 x 10 7 


0.17 


0.2 


15 


16 


1.04 


0.79 


0.81 


15.37 


15.6 


1.93 x 


10 6 


1.94 x 10 2 


3.20 x 


10 2 


0.49 


5.75 x 10 7 


0.16 


0.2 


15 


20 


1.08 


0.83 


0.86 


18.55 


16.2 


3.49 x 


10 6 


2.72 x 10 2 


4.26 x 


10 2 


0.48 


7.30 x 10 7 


0.17 


0.2 


20 


8 


0.98 


0.72 


0.68 


8.14 


19.7 


6.97 x 


10 5 


4.61 x 10 1 


9.66 x 


10 1 


0.39 


2.62 x 10 7 


0.23 


0.2 


20 


12 


1.05 


0.78 


0.76 


11.49 


20.9 


2.58 x 


10 6 


1.08 x 10 2 


1.99 x 


10 2 


0.49 


4.17 x 10 7 


0.18 


0.2 


20 


16 


1.07 


0.83 


0.81 


14.98 


21.4 


4.18 x 


10 6 


1.76 x 10 2 


3.06 x 


10 2 


0.49 


5.72 x 10 7 


0.16 


0.2 


20 


20 


1.10 


0.87 


0.86 


18.24 


21.9 


7.57 x 


10 6 


2.70 x 10 2 


4.57 x 


10 2 


0.52 


7.28 x 10 7 


0.16 


0.2 


25 


8 


0.98 


0.78 


0.67 


8.13 


24.6 


6.69 x 


10 5 


4.13 x 10 1 


1.06 x 


10 2 


0.25 


2.58 x 10 7 


0.25 


0.2 


25 


12 


1.05 


0.83 


0.77 


11.43 


26.3 


3.73 x 


10 6 


9.51 x 10 1 


1.82 x 


10 2 


0.45 


4.11 x 10 7 


0.18 


0.2 


25 


16 


1.09 


0.87 


0.82 


14.69 


27.2 


1.02 x 


10 7 


1.76 x 10 2 


3.10x 


10 2 


0.52 


5.64 x 10 7 


0.16 


0.2 


25 


20 


1.11 


0.90 


0.86 


18.08 


27.7 


1.59 x 


10 7 


2.53 x 10 2 


4.34 x 


10 2 


0.52 


7.18 x 10 7 


0.16 


0.3 


15 


8 


1.02 


0.68 


0.71 


7.83 


15.3 


1.43 x 


10" 


6.27 x 10 1 


1.02 x 


10 2 


0.64 


2.65 x 10 7 


0.26 


0.3 


15 


12 


1.04 


0.75 


0.77 


11.54 


15.6 


1.90 x 


10 6 


1.35 x 10 2 


2.19 x 


10 2 


0.56 


4.20 x 10 7 


0.21 


0.3 


15 


16 


1.08 


0.80 


0.83 


14.79 


16.2 


3.70 x 


10 6 


2.15 x 10 2 


3.44 x 


10 2 


0.58 


5.75 x 10 7 


0.19 


0.3 


15 


20 


1.13 


0.84 


0.89 


17.73 


16.9 


7.76 x 


10 6 


3.38 x 10 2 


5.42 x 


10 2 


0.61 


7.30 x 10 7 


0.20 


0.3 


20 


8 


1.02 


0.76 


0.70 


7.81 


20.5 


1.68 x 


10 6 


4.90 x 10 1 


8.94 x 


[()' 


0.45 


2.61 x 10 7 


0.29 


0.3 


20 


12 


1.08 


0.81 


0.79 


11.08 


21.7 


5.80 x 


10 6 


1.34 x 10 2 


2.23 x 


1 0- 


0.60 


4.14 x 10 7 


0.22 


0.3 


20 


16 


1.12 


0.84 


0.84 


14.29 


22.4 


1.25 x 


10 7 


2.20 x 10 2 


3.65 x 


1 0- 


0.63 


5.71 x 10 7 


0.19 


0.3 


20 


20 


1.14 


0.87 


0.88 


17.53 


22.8 


1.94 x 


10 7 


3.20 x 10 2 


5.26 x 


1 n2 


0.63 


7.27 x 10 7 


0.19 


0.3 


25 


8 


1.04 


0.82 


0.72 


7.66 


26.1 


3.17 x 


10 6 


5.10 x 10 1 


9.93 x 


[Q] 


0.42 


2.54 x 10 7 


0.30 


0.3 


25 


12 


1.09 


0.85 


0.80 


11.03 


27.2 


9.97 x 


10 6 


1.15 x 10 2 


1.91 x 


10 2 


0.55 


4.05 x 10 7 


0.23 


0.3 


25 


16 


1.13 


0.88 


0.85 


14.18 


28.2 


2.82 x 


10 7 


2.19 x 10 2 


3.49 x 


10 2 


0.63 


5.55 x 10 7 


0.20 


0.3 


25 


20 


1.16 


0.91 


0.90 


17.22 


29.1 


6.71 x 


10 7 


3.76 x 10 2 


5.84 x 


10 2 


0.68 


7.08 x 10 7 


0.21 


0.4 


15 


8 


1.04 


0.68 


0.72 


7.67 


15.6 


1.99 x 


10 6 


5.32 x 10 1 


8.93 x 


10 1 


0.55 


2.66 x 10 7 


0.32 


0.4 


15 


12 


1.07 


0.77 


0.80 


11.20 


16.1 


3.18x 


10" 


1.59 x 10 2 


2.50 x 


10 2 


0.63 


4.19 x 10 7 


0.24 


0.4 


15 


16 


1.11 


0.81 


0.85 


14.39 


16.7 


6.03 x 


10 6 


2.43 x 10 2 


3.89 x 


10 2 


0.63 


5.74 x 10 7 


0.21 


0.4 


15 


20 


1.16 


0.84 


0.90 


17.23 


17.4 


1.31 x 


10 7 


3.63 x 10 2 


5.79 x 


10 2 


0.67 


7.31 x 10 7 


0.22 


0.4 


20 


8 


1.06 


0.79 


0.74 


7.56 


21.2 


3.43 x 


10 6 


5.36 x 10 1 


7.34 x 


10 1 


0.52 


2.57 x 10 7 


0.34 


0.4 


20 


12 


1.11 


0.83 


0.82 


10.80 


22.2 


1.04 x 


10 7 


1.56 x 10 2 


2.38 x 


10 2 


0.64 


4.13 x 10 7 


0.25 


0.4 


20 


16 


1.16 


0.86 


0.89 


13.81 


23.2 


2.81 x 


10 7 


3.03 x 10 2 


4.55 x 


10 2 


0.72 


5.71 x 10 7 


0.24 


0.4 


20 


20 


1.20 


0.89 


0.93 


16.75 


23.9 


5.92 x 


10 7 


4.43 x 10 2 


6.68 x 


10 2 


0.73 


7.25 x 10 7 


0.25 


0.4 


25 


8 


1.09 


0.84 


0.75 


7.37 


27.1 


9.21 x 


10 6 


5.96 x 10 1 


1.00 x 


10 2 


0.45 


2.52 x 10 7 


0.36 


0.4 


25 


12 


1.11 


0.87 


0.84 


10.81 


27.8 


1.75 x 


10 7 


1.67 x 10 2 


2.26 x 


10 2 


0.61 


3.92 x 10 7 


0.27 


0.4 


25 


16 


1.16 


0.90 


0.90 


13.76 


29.1 


6.93 x 


10 7 


3.08 x 10 2 


4.16 x 


10 2 


0.71 


5.44 x 10 7 


0.25 


0.4 


25 


20 


1.21 


0.92 


0.95 


16.61 


30.1 


2.05 x 


10 8 


4.94 x 10 2 


6.45 x 


10 2 


0.72 


6.88 x 10 7 


0.28 
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Table A2. Numerical results for reference runs (7 > 0.5) 



7 


°0 


H* 


T* 

-■max 




IT'\ 


Nu 


a 

V 




v rms 


v * 


MM 






0.6 


10 


8 


0.95 


0.64 


0.70 


8.39 


9.5 


6.05 x 10 s 


9.53 x 


10 1 


1.36 x 10 2 


0.51 


2.68 x 10 7 


0.34 


0.6 


10 


12 


1.04 


0.71 


0.79 


11.55 


10.4 


1.54 x 10 6 


1.66 x 


LO 2 


2.41 x 10 2 


0.56 


4.23 x 10 7 


0.29 


0.6 


10 


16 


1.10 


0.76 


0.86 


14.52 


11.0 


3.07 x 10 6 


2.42 x 


10 2 


3.56 x 10 2 


0.56 


5.78 x 10 7 


0.28 


0.6 


10 


20 


1.16 


0.80 


0.92 


17.31 


11.6 


5.50 x 10 6 


3.15 x 


LO 2 


4.63 x 10 2 


0.55 


7.34 x 10 7 


0.28 


0.6 


15 


8 


1.10 


0.74 


0.76 


7.29 


16.5 


4.76 x 10 6 


6.65 x 


O 1 


1.03 x 10 2 


0.36 


2.61 x 10 7 


0.45 


0.6 


15 


12 


1.13 


0.79 


0.85 


10.60 


17.0 


8.32 x 10 6 


1.92 x 


O 2 


2.90 x 10 2 


0.69 


4.17 x 10 7 


0.30 


0.6 


15 


16 


1.18 


0.82 


0.90 


13.60 


17.7 


1.69 x 10 7 


2.84 x 


O 2 


4.39 x 10 2 


0.71 


5.71 x 10 7 


0.26 


0.6 


15 


20 


1.22 


0.85 


0.95 


16.35 


18.4 


3.52 x 10 7 


3.93 x 


10 2 


5.96 x 10 2 


0.71 


7.31 x 10 7 


0.27 


0.6 


20 


8 


1.13 


0.83 


0.79 


7.10 


22.5 


1.42 x 10 7 


6.93 x 


o 1 


8.92 x 10 1 


0.50 


2.54 x 10 7 


0.43 


0.6 


20 


12 


1.40 


1.38 


1.28 


8.59 


27.9 


3.91 x 10 9 


5.16x 


o 3 


9.91 x 10~ 2 


0.00 


3.41 x 10 7 


0.35 


0.6 


20 


16 


1.21 


0.88 


0.94 


13.20 


24.3 


8.69 x 10 7 


3.52 x 


LO 2 


4.85 x 10 2 


0.68 


5.63 x 10 7 


0.31 


0.6 


20 


20 


1.26 


0.89 


0.98 


15.85 


25.3 


2.45 x 10 s 


5.01 x 


o 2 


7.14 x 10 2 


0.66 


7.19 x 10 7 


0.32 


0.8 


10 


8 


1.09 


0.64 


0.79 


7.34 


10.9 


2.66 x 10 6 


1.05 x 


o 2 


1.57 x 10 2 


0.73 


2.69 x 10 7 


0.33 


0.8 


10 


12 


1.07 


0.72 


0.82 


11.17 


10.7 


2.28 x 10 6 


1.70 x 


10 2 


2.48 x 10 2 


0.58 


4.23 x 10 7 


0.32 


0.8 


10 


16 


1.14 


0.76 


0.88 


14.06 


11.4 


4.53 x 10 6 


2.44 x 


LO 2 


3.64 x 10 2 


0.58 


5.78 x 10 7 


0.30 


0.8 


10 


20 


1.20 


0.80 


0.94 


16.66 


12.0 


8.97 x 10 6 


3.30 x 


o 2 


5.00 x 10 2 


0.59 


7.36 x 10 7 


0.29 


0.8 


15 


8 


1.14 


0.77 


0.79 


7.03 


17.1 


8.92 x 10 6 


7.00 x 


o 1 


9.55 x 10 1 


0.33 


2.58 x 10 7 


0.50 


0.8 


15 


12 


1.18 


0.81 


0.89 


10.22 


17.6 


1.63 x 10 7 


1.93 x 


LO 2 


2.89 x 10 2 


0.64 


4.13 x 10 7 


0.32 


0.8 


15 


16 


1.23 


0.84 


0.95 


13.01 


18.5 


3.90 x 10 7 


2.98 x 


10 2 


4.43 x 10 2 


0.66 


5.74 x 10 7 


0.31 


0.8 


15 


20 


1.28 


0.86 


1.00 


15.59 


19.3 


9.11 x 10 7 


4.17 x 


o 2 


6.27 x 10 2 


0.66 


7.32 x 10 7 


0.32 


0.8 


20 


8 


1.32 


1.30 


1.18 


6.04 


26.5 


8.79 x 10 s 


1.79 x 


o 3 


4.93 x 10~ 3 


0.00 


2.06 x 10 7 


0.42 


0.8 


20 


12 


1.24 


0.88 


0.94 


9.69 


24.8 


1.47 x 10 8 


2.59 x 


0- 


2.90 x 10 2 


0.50 


4.04 x 10 7 


0.36 


0.8 


20 


16 


1.29 


0.89 


1.00 


12.38 


25.9 


4.51 x 10 8 


4.25 x 


o 2 


5.07 x 10 2 


0.53 


5.53 x 10 7 


0.34 


0.8 


20 


20 


1.35 


0.90 


1.05 


14.84 


27.0 


1.42 x 10 9 


6.75 x 


q2 


7.17 x 10 2 


0.52 


7.15 x 10 7 


0.31 


1.0 


10 


8 


1.10 


0.64 


0.79 


7.30 


11.0 


2.85 x 10 6 


1.00 x 


0- 


1.52 x 10 2 


0.60 


2.69 x 10 7 


0.39 


1.0 


10 


12 


1.11 


0.72 


0.84 


10.84 


11.1 


3.29 x 10 6 


1.67 x 


0^ 


2.50 x 10 2 


0.58 


4.24 x 10 7 


0.34 


1.0 


10 


16 


1.17 


0.76 


0.91 


13.64 


11.7 


6.67 x 10 6 


2.48 x 


q2 


3.80 x 10 2 


0.60 


5.80 x 10 7 


0.31 


1.0 


10 


20 


1.23 


0.79 


0.96 


16.21 


12.3 


1.29 x 10 7 


3.21 x 


0- 


5.02 x 10 2 


0.60 


7.38 x 10 7 


0.30 


1.0 


15 


8 


1.19 


0.86 


0.84 


6.74 


17.8 


1.94 x 10 7 


1.13 x 


o 2 


1.32 x 10 2 


0.27 


2.52 x 10 7 


0.55 


1.0 


15 


12 


1.24 


0.83 


0.91 


9.65 


18.7 


4.80 x 10 7 


1.86 x 


o 2 


2.79 x 10 2 


0.43 


4.14 x 10 7 


0.40 


1.0 


15 


16 


1.51 


1.49 


1.40 


10.61 


22.6 


3.08 x 10 9 


6.13x 


o 3 


4.93 x 10" 1 


0.00 


4.89 x 10 7 


0.30 


1.0 


15 


20 


1.52 


1.50 


1.42 


13.18 


22.8 


3.58 x 10 9 


8.31 x 


o 3 


1.55 x 10-' 


0.00 


6.32 x 10 7 


0.29 


1.0 


20 


8 


1.33 


1.31 


1.19 


6.04 


26.5 


8.89 x 10 8 


1.89 x 


o 3 


6.51 x 10~ 3 


0.00 


2.07 x 10 7 


0.43 


1.0 


20 


12 


1.39 


1.38 


1.28 


8.61 


27.9 


3.72 x 10 9 


4.71 x 


o 3 


4.37 x 10~ 3 


0.00 


3.41 x 10 7 


0.34 


1.0 


20 


16 


1.41 


1.40 


1.31 


11.38 


28.1 


4.76 x 10 9 


7.85 x 


o 3 


7.17 x 10- 3 


0.00 


4.82 x 10 7 


0.31 


1.0 


20 


20 


1.46 


1.45 


1.38 


13.69 


29.2 


1.49 x 10 10 


1.57 x 


I0 4 


2.07 x ur 3 


0.00 


6.30 x 10 7 


0.24 
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Table A3. Numerical results with shallow stiffening (7= 0.4, 9 = 15, and H* > 7) 



tfo 


H* 


h* 


A?) 


-"max 


IT*\ 
V 1 


Kin 


a 
u 


i\ui 




v rms 






Pn , 


*4' 






15 


8 


0.1 


3 


1.08 


0.73 


7.43 


16.1 


3.38 x 


10 6 


6.49 x 


10 1 


1.10 x 


10 2 


0.59 


2.64 x 


10 7 


0.38 


15 


8 


0.1 


10 


1.06 


0.75 


7.53 


15.9 


2.74 x 


10 6 


6.67 x 


10 1 


9.24 x 


10' 


0.71 


2.62 x 


10 7 


0.33 


15 


8 


0.1 


30 


1.05 


0.75 


7.62 


15.8 


2.28 x 


10 6 


7.51 x 


10 1 


1.10 x 


10 2 


0.69 


2.62 x 


10 7 


0.33 


15 


8 


0.1 


100 


1.07 


0.75 


7.49 


16.0 


3.01 x 


10 6 


7.62 x 


()' 


1.13 x 


10 2 


0.67 


2.62 x 


10 7 


0.33 


15 


8 


0.2 


3 


1.05 


0.76 


7.63 


15.7 


2.22 x 


10 6 


9.88 x 


o 1 


1.41 x 


10 2 


0.65 


2.60 x 


10 7 


0.33 


15 


8 


0.2 


10 


1.08 


0.79 


7.40 


16.2 


3.69 x 


10 6 


1.07 x 


o 2 


1.43 x 


10 2 


0.71 


2.59 x 


10 7 


0.33 


15 


8 


0.2 


30 


1.11 


0.82 


7.23 


16.6 


5.57 x 


10 6 


1.11 X 


10 2 


1.39 x 


10 2 


0.78 


2.58 x 


10 7 


0.31 


15 


8 


0.2 


100 


1.15 


0.87 


6.94 


17.3 


1.15 x 


10 7 


1.15 x 


10 2 


1.31 x 


10 2 


0.84 


2.57 x 


10 7 


0.28 


15 


8 


0.3 


3 


1.07 


0.77 


7.47 


16.1 


3.18x 


10 6 


9.47 x 


10 1 


1.34 x 


10 2 


0.62 


2.59 x 


10 7 


0.34 


15 


8 


0.3 


10 


1.11 


0.81 


7.21 


16.7 


5.80 x 


10 6 


1.05 x 


10 2 


1.34 x 


10 2 


0.65 


2.59 x 


10 7 


0.34 


15 


8 


0.3 


30 


1.15 


0.85 


6.98 


17.2 


1.03 x 


10 7 


1.07 x 


10 2 


1.27 x 


10 2 


0.65 


2.58 x 


10 7 


0.34 


15 


12 


0.1 


3 


1.08 


0.80 


11.13 


16.2 


3.50 x 


10 6 


1.65 x 


10 2 


2.55 x 


10 2 


0.65 


4.18 x 


10 7 


0.25 


15 


12 


0.1 


10 


1.10 


0.82 


10.94 


16.5 


4.77 x 


10 6 


1.89 x 


10 2 


2.87 x 


10 2 


0.69 


4.16 x 


10 7 


0.26 


15 


12 


0.1 


30 


1.10 


0.83 


10.87 


16.6 


5.29 x 


10 6 


1.87 x 


10 2 


2.79 x 


10 2 


0.71 


4.16 x 


10 7 


0.26 


15 


12 


0.1 


100 


1.12 


0.85 


10.74 


16.8 


6.54 x 


10 6 


1.89 x 


10 2 


2.73 x 


10 2 


0.74 


4.14 x 


10 7 


0.27 


15 


12 


0.2 


3 


1.11 


0.83 


10.86 


16.6 


5.42 x 


10 6 


1.91 x 


10 2 


2.82 x 


10 2 


0.67 


4.16x 


10 7 


0.26 


15 


12 


0.2 


10 


1.15 


0.88 


10.45 


17.2 


1.07 x 


10 7 


1.96 x 


o 2 


2.69 x 


10 2 


0.71 


4.16 x 


10 7 


0.27 


15 


12 


0.2 


30 


1.58 


1.43 


7.58 


23.7 


9.82 x 


10 9 


8.13x 


o 3 


1.57 x 10~ 2 


0.00 


3.09 x 


10 7 


0.36 


15 


12 


0.3 


3 


1.12 


0.84 


10.70 


16.8 


7.07 x 


10 6 


1.82 x 


o 2 


2.62 x 


10 2 


0.68 


4.17 x 


10 7 


0.26 


15 


12 


0.3 


10 


1.18 


0.90 


10.17 


17.7 


1.78 x 


10 7 


1.89 x 


10 2 


2.49 x 


10 2 


0.71 


4.16 x 


10 7 


0.27 


15 


12 


0.3 


30 


1.21 


0.93 


9.89 


18.2 


3.03 x 


10 7 


1.88 x 


10 2 


2.32 x 


10 2 


0.71 


4.15 x 


10 7 


0.29 


15 


16 


0.1 


3 


1.13 


0.87 


14.11 


17.0 


8.58 x 


10 6 


2.76 x 


10 2 


4.31 x 


10 2 


0.68 


5.74 x 


10 7 


0.23 


15 


16 


0.1 


10 


1.15 


0.89 


13.86 


17.3 


1.18 x 


10 7 


2.92 x 


o 2 


4.43 x 


10 2 


0.72 


5.74 x 


10 7 


0.24 


15 


16 


0.1 


30 


1.17 


0.91 


13.63 


17.6 


1.62 x 


10 7 


2.92 x 


o 2 


4.24 x 


10 2 


0.75 


5.72 x 


10 7 


0.25 


15 


16 


0.1 


100 


1.21 


0.95 


13.28 


18.1 


2.70 x 


10 7 


4.25 x 


o 2 


3.76 x 


10 2 


0.76 


5.63 x 


10 7 


0.27 


15 


16 


0.2 


3 


1.16 


0.89 


13.77 


17.4 


1.34 x 


10 7 


2.98 x 


10 2 


4.41 x 


10 2 


0.69 


5.74 x 


10 7 


0.24 


15 


16 


0.2 


10 


1.21 


0.94 


13.22 


18.2 


2.87 x 


10 7 


3.02 x 


10 2 


4.16 x 


10 2 


0.71 


5.73 x 


10 7 


0.26 


15 


16 


0.2 


30 


1.25 


0.97 


12.83 


18.7 


5.15 x 


10 7 


2.97 x 


10 2 


3.82 x 


10 2 


0.73 


5.71 x 


10 7 


0.28 


15 


16 


0.2 


100 


1.30 


1.02 


12.28 


19.6 


1.25 x 


10 8 


3.11 x 


10 2 


3.80 x 


10 2 


0.72 


5.71 x 


10 7 


0.29 


15 


16 


0.3 


3 


1.18 


0.91 


13.53 


17.8 


1.86 x 


10' 


2.95 x 


10 2 


4.23 x 


10 2 


0.72 


5.73 x 


10' 


0.24 


15 


16 


0.3 


10 


1.24 


0.96 


12.96 


18.5 


4.21 x 


10 7 


2.81 x 


10 2 


3.64 x 


10 2 


0.74 


5.73 x 


10 7 


0.27 


15 


16 


0.3 


30 


1.28 


0.99 


12.47 


19.3 


9.10 x 


10 7 


2.84 x 


o 2 


3.46 x 


10 2 


0.75 


5.72 x 


10 7 


0.28 


15 


16 


0.3 


100 


1.33 


1.03 


12.04 


20.0 


1.89 x 


10 s 


2.86 x 


o 2 


3.25 x 


10 2 


0.74 


5.69 x 


10 7 


0.29 


15 


20 


0.1 


3 


1.18 


0.93 


16.89 


17.8 


1.89 x 


10 7 


3.89 x 


o 2 


6.00 x 


10 2 


0.70 


7.31 x 


10 7 


0.23 


15 


20 


0.1 


10 


1.21 


0.95 


16.55 


18.1 


2.79 x 


10 7 


3.87 x 


10 2 


5.70 x 


10 2 


0.73 


7.30 x 


10 7 


0.25 


15 


20 


0.1 


100 


1.27 


1.00 


15.81 


19.0 


6.82 x 


10 7 


3.94 x 


10 2 


5.42 x 


10 2 


0.78 


7.29 x 


10 7 


0.28 


15 


20 


0.2 


3 


1.21 


0.95 


16.53 


18.2 


2.85 x 


10 7 


3.93 x 


10 2 


5.73 x 


10 2 


0.70 


7.31 x 


10 7 


0.25 


15 


20 


0.2 


10 


1.25 


0.99 


15.96 


18.8 


5.70 x 


10 7 


3.99 x 


10 2 


5.35 x 


10 2 


0.71 


7.30 x 


10 7 


0.27 


15 


20 


0.2 


30 


1.30 


1.02 


15.45 


19.4 


1.08 x 


10 s 


3.99 x 


10 2 


5.12 x 


10 2 


0.72 


7.29 x 


10 7 


0.28 


15 


20 


0.3 


3 


1.22 


0.96 


16.40 


18.3 


3.34 x 


10 7 


3.80 x 


10 2 


5.36 x 


10 2 


0.71 


7.30 x 


10 7 


0.25 


15 


20 


0.3 


10 


1.28 


1.01 


15.67 


19.2 


8.20 x 


10 7 


3.78 x 


10 2 


4.85 x 


10 2 


0.73 


7.29 x 


10 7 


0.28 


15 


20 


0.3 


30 


1.33 


1.05 


15.07 


19.9 


1.83 x 


10 s 


3.94 x 


10 2 


4.69 x 


10 2 


0.74 


7.27 x 


10 7 


0.30 


15 


20 


0.3 


100 


1.39 


1.09 


14.36 


20.9 


5.13 x 


10 8 


4.03 x 


10 2 


4.53 x 


10 2 


0.76 


7.25 x 


10 7 


0.30 
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Table A4. Numerical results with shallow stiffening (y = 0.4, 9 = 20, and H* > 7) 



fo 


n 


ft 


Atj 


'max 







u 




v rms 


v., 




r()A 


VP 






20 


8 


0.1 


3 


1.06 


0.74 


7.52 


21.3 


3.80 x 10 6 


5.14 x 


10 1 


7.50 x 


10 1 


0.51 


2.58 x 


10 7 


0.34 


20 


8 


0.1 


10 


1.08 


0.74 


7.42 


21.6 


5.10 x 10 6 


5.77 x 


O 1 


9.57 x 


10' 


0.50 


2.57 x 


10 7 


0.37 


20 


8 


0.1 


30 


1.09 


0.75 


7.32 


21.9 


7.02 x 10 6 


5.60 x 


o 1 


9.39 x 


10 1 


0.49 


2.56 x 


10 7 


0.38 


20 


8 


0.1 


100 


1.12 


0.76 


7.11 


22.5 


1.37 x 10 7 


5.59 x 


o 1 


8.57 x 


10 1 


0.52 


2.52 x 


10 7 


0.38 


20 


8 


0.2 


3 


1.10 


0.75 


7.28 


22.0 


7.94 x 10 6 


5.74 x 


10 1 


9.07 x 


10' 


0.47 


2.57 x 


10 7 


0.37 


20 


8 


0.2 


10 


1.13 


0.78 


7.08 


22.6 


1.52 x 10 7 


5.34 x 


10 1 


7.60 x 


10' 


0.50 


2.57 x 


10 7 


0.34 


20 


8 


0.2 


30 


1.14 


0.81 


7.02 


22.8 


1.89 x 10 7 


7.63 x 


o 1 


9.28 x 


10' 


0.63 


2.53 x 


10 7 


0.33 


20 


8 


0.3 


3 


1.10 


0.76 


7.29 


22.0 


7.81 x 10 6 


5.88 x 


o 1 


9.37 x 


10' 


0.48 


2.57 x 


10 7 


0.37 


20 


8 


0.3 


10 


1.13 


0.80 


7.06 


22.7 


1.62 x 10 7 


6.88 x 


10 1 


9.41 x 


10' 


0.53 


2.54 x 


10 7 


0.36 


20 


8 


0.3 


30 


1.15 


0.85 


6.95 


23.0 


2.42 x 10 7 


1.17x 


o 2 


1.30 x 


10 2 


0.62 


2.48 x 


10 7 


0.36 


20 


8 


0.3 


100 


1.20 


0.90 


6.69 


23.9 


6.19 x 10 7 


1.32 x 


o 2 


1.24 x 


10 2 


0.65 


2.45 x 


10 7 


0.34 


20 


12 


0.1 


3 


1.12 


0.82 


10.74 


22.4 


1.18 x 10 7 


1.50 x 


o 2 


2.32 x 


10 2 


0.65 


4.12 x 


10 7 


0.26 


20 


12 


0.1 


10 


1.13 


0.84 


10.65 


22.5 


1.44 x 10 7 


1.77 x 


o 2 


2.58 x 


10 2 


0.69 


4.08 x 


10 7 


0.27 


20 


12 


0.1 


30 


1.14 


0.86 


10.56 


22.7 


1.75 x 10 7 


2.01 x 


10 2 


2.72 x 


10 2 


0.72 


4.02 x 


10 7 


0.28 


20 


12 


0.1 


100 


1.15 


0.89 


10.42 


23.1 


2.47 x 10 7 


2.03 x 


o 2 


2.48 x 


10 2 


0.74 


3.97 x 


10 7 


0.29 


20 


12 


0.2 


3 


1.14 


0.86 


10.54 


22.8 


1.85 x 10 7 


2.05 x 


o 2 


2.86 x 


10 2 


0.70 


4.09 x 


10 7 


0.27 


20 


12 


0.2 


10 


1.17 


0.89 


10.30 


23.3 


3.22 x 10 7 


2.23 x 


10 2 


2.86 x 


10 2 


0.72 


4.05 x 


10 7 


0.29 


20 


12 


0.2 


30 


1.19 


0.92 


10.08 


23.8 


5.51 x 10 7 


2.33 x 


o 2 


2.68 x 


10 2 


0.75 


4.03 x 


10 7 


0.32 


20 


12 


0.3 


3 


1.15 


0.86 


10.43 


23.0 


2.36 x 10 7 


1.96 x 


10 2 


2.74 x 


10 2 


0.69 


4.09 x 


10 7 


0.28 


20 


12 


0.3 


10 


1.18 


0.90 


10.18 


23.6 


4.26 x 10 7 


2.05 x 


10 2 


2.57 x 


10 2 


0.71 


4.06 x 


10 7 


0.29 


20 


12 


0.3 


30 


1.21 


0.93 


9.91 


24.2 


8.38 x 10 7 


2.13 x 


o 2 


2.46 x 


10 2 


0.71 


4.03 x 


10 7 


0.31 


20 


12 


0.3 


100 


1.25 


0.96 


9.59 


25.0 


1.94 x 10 8 


2.19 x 


10 2 


2.30 x 


10 2 


0.71 


3.94 x 


10 7 


0.33 


20 


16 


0.1 


3 


1.17 


0.89 


13.73 


23.3 


3.22 x 10 7 


3.15 x 


o 2 


4.66 x 


10 2 


0.73 


5.66 x 


10 7 


0.25 


20 


16 


0.1 


10 


1.17 


0.91 


13.64 


23.5 


3.79 x 10 7 


3.22 x 


o 2 


4.48 x 


10 2 


0.75 


5.61 x 


10 7 


0.26 


20 


16 


0.1 


30 


1.18 


0.92 


13.54 


23.6 


4.52 x 10 7 


3.27 x 


o 2 


4.28 x 


10 2 


0.75 


5.55 x 


10 7 


0.28 


20 


16 


0.1 


100 


1.21 


0.95 


13.19 


24.3 


8.75 x 10 7 


3.38 x 


o 2 


4.03 x 


10 2 


0.78 


5.50 x 


10 7 


0.29 


20 


16 


0.2 


3 


1.18 


0.91 


13.52 


23.7 


4.69 x 10 7 


3.31 x 


o 2 


4.62 x 


10 2 


0.72 


5.66 x 


10 7 


0.26 


20 


16 


0.2 


10 


1.22 


0.95 


13.18 


24.3 


8.98 x 10 7 


3.47 x 


o 2 


4.38 x 


10 2 


0.74 


5.61 x 


10 7 


0.29 


20 


16 


0.2 


30 


1.25 


0.98 


12.81 


25.0 


1.87 x 10 8 


3.64 x 


o 2 


4.34 x 


1() 2 


0.74 


5.58 x 


10 7 


0.31 


20 


16 


0.2 


100 


1.29 


1.02 


12.41 


25.8 


4.24 x 10 8 


3.82 x 


o 2 


3.94 x 


10 2 


0.73 


5.55 x 


10 7 


0.34 


20 


16 


0.3 


3 


1.19 


0.92 


13.46 


23.8 


5.27 x 10 7 


3.15 x 


o 2 


4.38 x 


w l 


0.72 


5.64 x 


10' 


0.26 


20 


16 


0.3 


10 


1.23 


0.96 


13.00 


24.6 


1.26 x 10 8 


3.23 x 


10 2 


4.01 x 


10 2 


0.74 


5.60 x 


10 7 


0.28 


20 


16 


0.3 


30 


1.27 


0.99 


12.60 


25.4 


2.89 x 10 8 


3.32 x 


10 2 


3.80 x 


10 2 


0.75 


5.58 x 


10 7 


0.30 


20 


16 


0.3 


100 


1.31 


1.02 


12.23 


26.2 


6.45 x 10 8 


3.34 x 


o 2 


3.47 x 


10 2 


0.75 


5.55 x 


10 7 


0.32 


20 


20 


0.1 


3 


1.21 


0.95 


16.54 


24.2 


8.04 x 10 7 


4.53 x 


10 2 


6.49 x 


10 2 


0.75 


7.23 x 


10 7 


0.27 


20 


20 


0.1 


10 


1.22 


0.96 


16.37 


24.5 


1.05 x 10 8 


4.65 x 


o 2 


6.32 x 


10 2 


0.77 


7.18 x 


10 7 


0.28 


20 


20 


0.1 


100 


1.29 


1.01 


15.48 


25.9 


4.59 x 10 8 


4.89 x 


o 2 


6.44 x 


10 2 


0.82 


7.06 x 


10 7 


0.27 


20 


20 


0.2 


3 


1.23 


0.96 


16.30 


24.6 


1.17 x 10 8 


4.66 x 


o 2 


6.42 x 


10 2 


0.74 


7.22 x 


10 7 


0.28 


20 


20 


0.2 


10 


1.26 


0.99 


15.88 


25.2 


2.31 x 10 8 


4.70 x 


o 2 


5.94 x 


10 2 


0.74 


7.19 x 


10 7 


0.30 


20 


20 


0.2 


30 


1.30 


1.02 


15.42 


26.0 


5.07 x 10 8 


4.90 x 


o 2 


5.88 x 


10 2 


0.74 


7.17 x 


10 7 


0.31 


20 


20 


0.2 


100 


1.35 


1.06 


14.88 


26.9 


1.34 x 10 9 


5.16x 


o 2 


5.64 x 


10 2 


0.74 


7.13 x 


10 7 


0.31 


20 


20 


0.3 


3 


1.23 


0.97 


16.25 


24.6 


1.28 x 10 8 


4.44 x 


10 2 


6.03 x 


10 2 


0.74 


7.22 x 


10 7 


0.28 


20 


20 


0.3 


10 


1.28 


1.00 


15.68 


25.5 


3.27 x 10 8 


4.45 x 10 2 


5.51 x 


10 2 


0.76 


7.18 x 


10 7 


0.29 


20 


20 


0.3 


100 


1.36 


1.06 


14.78 


27.1 


1.65 x 10 9 


4.59 x 


10 2 


4.78 x 


10 2 


0.76 


7.13 x 


10 7 


0.33 
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Table A5. Numerical results with shallow stiffening (y = 0.4, 9 = 25, and H* > 7) 



00 


ft 


n 


A?] 


-"max 




Nu 


n 
U 


Raj 


v rms 




Poa 


<P 


, 1 , ,fi 

1> S /9 


25 


8 


0.1 


3 


1.09 


0.76 


7.37 


21.1 


9.32 x 10 6 


6.88 x 10 1 


1.07 x 10 2 


0.49 


2.46 x 10 7 


0.36 


25 


8 


0.1 


10 


1.29 


1.16 


6.20 


32.3 


1.86 x 10 9 


2.99 x 10 3 


6.71 x 10~ 4 


0.00 


2.05 x 10 7 


0.43 


25 


8 


0.1 


30 


1.29 


1.16 


6.22 


32.2 


1.67 x 10 9 


3.02 x 10 3 


9.93 x 10~ 5 


0.00 


2.06 x 10 7 


0.43 


25 


8 


0.1 


100 


1.11 


0.80 


7.20 


27.8 


1.87 x 10 7 


8.04 x 10 1 


9.53 x 10 1 


0.55 


2.40 x 10 7 


0.34 


25 


8 


0.2 


3 


1.11 


0.77 


7.22 


27.7 


1.69 x 10 7 


6.59 x 10 1 


9.62 x 10 1 


0.49 


2.48 x 10 7 


0.36 


25 


8 


0.2 


10 


1.14 


0.80 


7.00 


28.6 


4.12 x 10 7 


7.04 x 10' 


9.08 x 10 1 


0.54 


2.45 x 10 7 


0.35 


25 


8 


0.2 


30 


1.20 


0.86 


6.66 


30.0 


1.84 x 10 8 


1.00 x 10 2 


6.78 x 10 1 


0.47 


2.47 x 10 7 


0.32 


25 


8 


0.2 


100 


1.20 


0.88 


6.70 


29.9 


1.60 x 10* 


1.62 x 10 2 


8.69 x 10 1 


0.66 


2.35 x 10 7 


0.32 


25 


8 


0.3 


3 


1.10 


0.78 


7.25 


27.6 


1.46 x 10 7 


6.60 x 10 1 


9.81 x 10 1 


0.49 


2.46 x 10 7 


0.36 


25 


12 


0.1 


3 


1.13 


0.84 


10.60 


28.3 


3.10 x 10 7 


1.42 x 10 2 


2.11 x 10 2 


0.62 


3.93 x 10 7 


0.28 


25 


12 


0.1 


10 


1.14 


0.85 


10.50 


28.6 


4.14 x 10 7 


1.60 x 10 2 


2.19 x 10 2 


0.66 


3.88 x 10 7 


0.28 


25 


12 


0.1 


30 


1.16 


0.89 


10.34 


29.0 


6.60 x 10 7 


2.09 x 10 2 


2.34 x 10 2 


0.72 


3.80 x 10 7 


0.28 


25 


12 


0.3 


3 


1.15 


0.86 


10.42 


28.8 


5.22 x 10 7 


1.72 x 10 2 


2.28 x 10 2 


0.66 


3.90 x 10 7 


0.28 


25 


12 


0.3 


10 


1.18 


0.90 


10.14 


29.6 


1.20 x 10 s 


2.12 x 10 2 


2.42 x 10 2 


0.70 


3.79 x 10 7 


0.31 


25 


16 


0.1 


30 


1.22 


0.95 


13.09 


30.6 


3.33 x 10 8 


3.71 x 10 2 


4.17 x 10 2 


0.76 


5.42 x 10 7 


0.28 


25 


16 


0.2 


3 


1.19 


0.93 


13.41 


29.8 


1.52 x 10 8 


3.62 x 10 2 


4.41 x 10 2 


0.74 


5.32 x 10 7 


0.28 


25 


16 


0.2 


10 


1.22 


0.96 


13.12 


30.5 


3.08 x 10 8 


3.86 x 10 2 


4.11 x 10 2 


0.74 


5.72 x 10 7 


0.29 


25 


16 


0.3 


10 


1.23 


0.96 


13.04 


30.7 


3.72 x 10 8 


6.10 x 10 2 


4.00 x 10 2 


0.76 


7.53 x 10 7 


0.37 


25 


20 


0.1 


3 


1.22 


0.96 


16.35 


30.6 


3.38 x 10 8 


5.04 x 10 2 


6.82 x 10 2 


0.76 


6.81 x 10 7 


0.27 


25 


20 


0.1 


10 


1.25 


0.97 


15.96 


31.4 


7.30 x 10 8 


5.22 x 10 2 


7.03 x 10 2 


0.78 


6.79 x 10 7 


0.27 


25 


20 


0.1 


30 


1.30 


0.97 


15.33 


32.6 


2.64 x 10 9 


4.88 x 10 2 


7.25 x 10 2 


0.77 


7.06 x 10 7 


0.22 


25 


20 


0.1 


100 


1.35 


1.06 


14.77 


33.9 


9.60 x 10 9 


1.28 x 10 3 


6.47 x 10 2 


0.80 


6.99 x 10 7 


0.21 
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Table A6. Numerical results with shallow stiffening (7= 0.4, 9 = 15, and H* < 7) 





H* 


/,* 


Atj 


T* 

-"max 


/T*\ 
\* I 


Nil 


a 
u 




•rms 




D.1 






15 


2 


0.1 


3 


0.82 


0.51 


2.44 


12.3 


5.59 x 10 4 


7.67 x 10° 


5.19 x 10° 


0.15 


2.44 x 10 6 


0.88 


15 


2 


0.1 


10 


0.82 


0.51 


2.44 


12.3 


5.59 x 10 4 


7.67 x 10° 


5.18 x 10° 


0.20 


2.44 x 10 6 


0.88 


15 


2 


0.1 


30 


0.82 


0.51 


2.44 


12.3 


5.59 x 10 4 


7.67 x 10° 


5.17 x 10° 


0.25 


2.44 x 10 6 


0.88 


15 


2 


0.1 


100 


0.82 


0.51 


2.44 


12.3 


5.59 x 10 4 


7.67 x 10° 


5.19 x 10° 


0.27 


2.44 x 10 6 


0.88 


15 


2 


0.2 


3 


0.83 


0.52 


2.41 


12.4 


6.44 x 10 4 


7.88 x 10° 


4.85 x 10° 


0.35 


2.42 x 10 6 


0.88 


15 


2 


0.2 


10 


0.84 


0.52 


2.38 


12.6 


7.50 x 10 4 


8.17 x 10° 


4.80 x 10° 


0.35 


2.54 x 10 6 


0.88 


15 


2 


0.2 


30 


0.84 


0.52 


2.38 


12.6 


7.70 x 10 4 


8.20 x 10° 


4.75 x 10° 


0.45 


2.54 x 10 6 


0.87 


15 


2 


0.2 


100 


0.84 


0.52 


2.38 


12.6 


7.74 x 10 4 


8.21 x 10° 


4.74 x 10° 


0.52 


2.55 x 10 6 


0.87 


15 


2 


0.3 


3 


0.85 


0.52 


2.36 


12.7 


8.78 x 10 4 


8.27 x 10° 


4.55 x 10° 


0.38 


2.50 x 10 6 


0.87 


15 


2 


0.3 


10 


0.86 


0.54 


2.31 


13.0 


1.13 x 10 5 


8.76 x 10° 


3.94 x 10° 


0.46 


2.36 x 10 6 


0.89 


15 


2 


0.3 


30 


0.91 


0.63 


2.19 


13.7 


2.54 x 10 5 


1.18 x 10 1 


4.22 x 10~ 3 


0.00 


1.17 x 10 6 


0.87 


15 


2 


0.3 


100 


0.91 


0.63 


2.19 


13.7 


2.55 x 10 s 


1.18 x 10 1 


1.24 x 10~ 3 


0.00 


1.16 x 10 6 


0.87 


15 


4 


0.1 


3 


0.90 


0.60 


4.42 


13.6 


2.16 x 10 5 


2.28 x 10 1 


2.60 x 10 1 


0.50 


1.07 x 10 7 


0.52 


15 


4 


0.1 


10 


0.91 


0.60 


4.42 


13.6 


2.19 x 10 5 


2.28 x 10 1 


2.59 x 10 1 


0.45 


1.07 x 10 7 


0.52 


15 


4 


0.1 


30 


0.91 


0.60 


4.40 


13.6 


2.34 x 10 s 


2.28 x 10 1 


2.59 x 10 1 


0.54 


1.07 x 10 7 


0.52 


15 


4 


0.1 


100 


0.91 


0.60 


4.40 


13.7 


2.36 x 10 s 


2.28 x 10 1 


2.59 x 10 1 


0.55 


1.07 x 10 7 


0.51 


15 


4 


0.2 


3 


0.91 


0.61 


4.38 


13.7 


2.49 x 10 s 


2.26 x 10 1 


2.50 x 10 1 


0.48 


1.07 x 10 7 


0.49 


15 


4 


0.2 


10 


0.95 


0.62 


4.23 


14.2 


4.22 x 10 5 


2.58 x 10 1 


2.98 x 10 1 


0.44 


1.05 x 10 7 


0.52 


15 


4 


0.2 


30 


0.94 


0.63 


4.23 


14.2 


4.12 x 10 5 


2.54 x 10 1 


2.53 x 10 1 


0.70 


1.05 x 10 7 


0.48 


15 


4 


0.2 


100 


0.95 


0.64 


4.19 


14.3 


4.85 x 10 5 


2.58 x 10 1 


2.45 x 10 1 


0.71 


1.05 x 10 7 


0.48 


15 


4 


0.3 


3 


0.92 


0.62 


4.33 


13.9 


2.95 x 10 s 


2.28 x 10 1 


2.35 x 10 1 


0.46 


1.06 x 10 7 


0.48 


15 


4 


0.3 


10 


0.94 


0.65 


4.24 


14.2 


4.09 x 10 5 


2.33 x 10 1 


2.12 x 10 1 


0.65 


1.05 x 10 7 


0.46 


15 


4 


0.3 


30 


0.96 


0.67 


4.17 


14.4 


5.25 x 10 5 


2.41 X 10 


1 nr\ * , 1 aI 

1.99 X 10 


0.55 


1 r\A . , ia7 

1.04 X 10 


0.46 


15 


4 


0.3 


100 


0.99 


0.69 


4.06 


14.8 


8.02 x 10 s 


2.62 x 10 1 


2.05 x 10 1 


0.54 


1.02 x 10 7 


0.49 


15 


6 


0.1 


3 


0.99 


0.67 


6.07 


14.8 


8.22 x 10 s 


3.67 x 10 1 


5.24 x 10 1 


0.62 


1.86 x 10 7 


0.39 


15 


6 


0.1 


10 


0.99 


0.68 


6.05 


14.9 


8.74 x 10 5 


3.65 x 10 1 


5.15 x 10 1 


0.65 


1.86 x 10 7 


0.38 


15 


6 


0.1 


30 


0.99 


0.68 


6.04 


14.9 


9.09 x 10 5 


3.64 x 10 1 


5.08 x 10 1 


0.67 


1.86 x 10 7 


0.38 


15 


6 


0.1 


100 


1.02 


0.68 


5.89 


15.3 


1.34 x 10 6 


3.91 x 10 1 


5.51 x 10 1 


0.60 


1.85 x 10 7 


0.40 


15 


6 


0.2 


3 


1.00 


0.69 


6.02 


15.0 


9.50 x 10 5 


3.66 x 10 1 


4.95 x 10 1 


0.66 


1.86 x 10 7 


0.36 


15 


6 


0.2 


10 


1.03 


0.71 


5.80 


15.5 


1.73 xlO 6 


4.14 x 10 1 


5.46 x 10 1 


0.71 


1.84 x 10 7 


0.35 


15 


6 


0.2 


30 


1.06 


0.74 


5.65 


15.9 


2.72 x 10 6 


4.60 x 10 1 


5.43 x 10 1 


0.75 


1.83 x 10 7 


0.33 


15 


6 


0.2 


100 


1.10 


0.77 


5.45 


16.5 


5.06 x 10 6 


5.91 x 10 1 


6.26 x 10 1 


0.82 


1.80 x 10 7 


0.35 


15 


6 


0.3 


3 


1.01 


0.70 


5.93 


15.2 


1.21 x 10 6 


3.65 x 10 1 


4.78 x 10 1 


0.66 


1.85 x 10 7 


0.37 


15 


6 


0.3 


10 


1.08 


0.76 


5.57 


16.2 


3.40 x 10 6 


4.68 x 10 1 


5.33 x 10 1 


0.59 


1.81 x 10 7 


0.39 


15 


6 


0.3 


30 


1.11 


0.80 


5.39 


16.7 


6.03 x 10 6 


7.22 x 10 1 


7.49 x 10 1 


0.63 


1.78 x 10 7 


0.40 


15 


6 


0.3 


100 


1.16 


0.84 


5.18 


17.4 


1.24 x 10 7 


7.88 x 10 1 


6.61 x 10 1 


0.68 


1.76 x 10 7 


0.36 
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Table A7. Numerical results with shallow stiffening (7 = 0.4, 9 = 20 & 25, and H* < 7) 





H* 


/,* 


Atj 


T* 

-"max 


/T*\ 
\* I 




a 
u 






v rms 




D.l 








20 


2 


0.1 


3 


0.96 


0.65 


2.08 


19.3 


4.62 x 


10 s 


9.88 x 


10° 


2.53 x 10~ 3 


0.00 


7.02 x 


10 5 


0.94 


20 


2 


0.2 


3 


0.96 


0.65 


2.08 


19.3 


4.62 x 


io 5 


9.89 x 


10° 


1.02 x 10~ 3 


0.00 


7.02 x 


IO 5 


0.94 


20 


2 


0.3 


3 


0.96 


0.65 


2.08 


19.3 


4.66 x 


io 5 


9.90 x 


10° 


8.63 x Kr 4 


0.00 


6.98 x 


IO 5 


0.94 


20 


4 


0.1 


3 


1.01 


0.66 


3.97 


20.2 


1.18 x 


IO 6 


2.28 x 


IO 1 


2.09 x IO 1 


0.56 


9.79 x 


IO 6 


0.58 


20 


4 


0.1 


10 


0.98 


0.64 


4.07 


19.7 


7.25 x 


IO 5 


2.21 x 


IO 1 


4.45 x IO 1 


0.37 


9.86 x 


IO 6 


0.53 


20 


4 


0.1 


30 


1.01 


0.67 


3.94 


20.3 


1.35 x 


IO 6 


2.33 x 


10' 


2.03 x IO 1 


0.58 


9.66 x 


IO 6 


0.57 


20 


4 


0.1 


100 


1.02 


0.67 


3.93 


20.4 


1.45 x 


IO 6 


2.35 x 


IO 1 


2.08 x IO 1 


0.50 


9.59 x 


IO 6 


0.56 


20 


4 


0.2 


3 


1.00 


0.66 


3.99 


20.1 


1.09 x 


IO 6 


2.62 x 


IO 1 


4.06 x IO 1 


0.43 


9.70 x 


IO 6 


0.52 


20 


4 


0.2 


10 


1.02 


0.67 


3.93 


20.4 


1.46 x 


IO 6 


2.57 x 


IO 1 


3.53 x IO 1 


0.53 


9.70 x 


IO 6 


0.50 


20 


4 


0.2 


30 


1.04 


0.69 


3.83 


20.9 


2.56 x 


IO 6 


2.47 x 


IO 1 


2.05 x IO 1 


0.49 


9.70 x 


IO 6 


0.54 


20 


4 


0.2 


100 


1.05 


0.69 


3.81 


21.0 


2.86 x 


IO 6 


2.67 x 


IO 1 


2.63 x IO 1 


0.50 


9.57 x 


IO 6 


0.53 


20 


4 


0.3 


3 


1.03 


0.69 


3.89 


20.6 


1.82 x 


IO 6 


2.38 x 


IO 1 


1.87 x IO 1 


0.58 


9.74 x 


IO 6 


0.5 1 


20 


4 


0.3 


10 


1.21 


0.99 


3.31 


24.1 


7.64 x 


IO 7 


3.69 x 


10 2 


9.90 x 10~ 7 


0.00 


6.30 x 


IO 6 


0.60 


20 


6 


0.1 


3 


1.05 


0.71 


5.73 


20.9 


2.66 x 


IO 6 


3.78 x 


IO 1 


4.56 x IO 1 


0.49 


1.78 x 


IO 7 


0.4 1 


20 


6 


0.1 


10 


1.04 


0.70 


5.78 


20.8 


2.28 x 


IO 6 


3.73 x 


IO 1 


5.86 x IO 1 


0.45 


1.78 x 


IO 7 


0.43 


20 


6 


0.1 


30 


1.07 


0.72 


5.63 


21.3 


3.99 x 


IO 6 


3.48 x 


IO 1 


4.54 x IO 1 


0.52 


1.77 x 


IO 7 


0.40 


20 


6 


0.1 


100 


1.08 


0.72 


5.58 


21.5 


4.85 x 


IO 6 


3.58 x 


IO 1 


5.00 x IO 1 


0.46 


1.76 x 


IO 7 


0.43 


20 


6 


0.2 


3 


1.06 


0.71 


5.65 


21.2 


3.66 x 


IO 6 


3.35 x 


IO 1 


4.62 x IO 1 


0.46 


1.78 x 


IO 7 


0.42 


20 


6 


0.2 


10 


1.09 


0.73 


5.51 


21.8 


6.55 x 


IO 6 


3.72 x 


IO 1 


4.73 x IO 1 


0.53 


1.77 x 


IO 7 


0.42 


20 


6 


0.2 


30 


1.09 


0.74 


5.51 


21.8 


6.36 x 


IO 6 


4.24 x 


10 


5.40 x IO 1 


0.62 


1.76 x 


IO 7 


0.41 


20 


6 


0.2 


100 


1.13 


0.78 


5.31 


22.6 


1.51 x 


10 


4.90 x 


IO 1 


4.60 X 10 


0.68 


1.74 x 


10 


0.37 


20 


6 


0.3 


3 


1.07 


0.72 


5.63 


21.3 


3.93 x 


IO 6 


3.78 x 


IO 1 


5.60 x IO 1 


0.42 


1.77 x 


IO 7 


0.46 


20 


6 


0.3 


10 


1.11 


0.76 


5.42 


22.1 


9.33 x 


IO 6 


3.96 x 


10' 


4.84 x IO 1 


0.45 


1.76 x 


IO 7 


0.42 


20 


6 


0.3 


30 


1.14 


0.79 


5.27 


22.8 


1.80 x 


IO 7 


4.60 x 


IO 1 


4.96 x IO 1 


0.57 


1.74 x 


IO 7 


0.37 


20 


6 


0.3 


100 


1.43 


1.21 


4.20 


28.6 


7.48 x 


IO 9 


3.65 x 


IO 3 


1.97 x 10~ 8 


0.00 


8.70 x 


IO 6 


0.53 


25 


2 


0.1 


3 


0.98 


0.66 


2.04 


24.6 


6.29 x 


IO 5 


7.35 x 


10° 


5.99 x 10~ 8 


0.00 


3.93 x 


IO 5 


0.96 


25 


2 


0.2 


3 


0.98 


0.66 


2.04 


24.6 


6.29 x 


IO 5 


7.35 x 


10° 


5.48 x 10~ 8 


0.00 


3.93 x 


IO 5 


0.96 


25 


2 


0.3 


3 


0.98 


0.66 


2.04 


24.6 


6.29 x 


IO 5 


7.36 x 


10° 


5.06 x IO" 8 


0.00 


3.93 x 


IO 5 


0.96 


25 


4 


0.1 


3 


1.17 


0.97 


3.41 


29.3 


8.56 x 


IO 7 


3.90 x 


10 2 


8.75 x 10~ 6 


0.00 


6.68 x 


IO 6 


0.60 


25 


4 


0.2 


3 


1.17 


0.97 


3.41 


29.3 


8.59 x 


IO 7 


3.90 x 


10 2 


4.69 x 10~ 6 


0.00 


6.66 x 


IO 6 


0.60 


25 


4 


0.3 


3 


1.19 


0.98 


3.36 


29.7 


1.33 x 


IO 8 


4.66 x 


10 2 


4.88 x 10~ 8 


0.00 


6.40 x 


IO 6 


0.60 


25 


6 


0.1 


3 


1.25 


1.09 


4.80 


31.2 


6.39 x 


IO 8 


l.39x 


IO 3 


1.79 x 10~ 4 


0.00 


1.36 x 


IO 7 


0.49 


25 


6 


0.2 


3 


1.25 


1.09 


4.79 


31.3 


7.05 x 


IO 8 


1.23 x 


IO 3 


2.44 x 10~ 5 


0.00 


1.33 x 


IO 7 


0.46 


25 


6 


0.3 


3 


1.28 


1.12 


4.67 


32.1 


1.56 x 


IO" 


2.00 x 


IO 3 


3.59 x IO" 4 


0.00 


1.31 x 


IO 7 


0.50 



51 



Table A8. Numerical results with shallow stiffening (y = 0.8) 



On 
"0 


H* 


h* 


Arj 




!T*\ 


Nu 


Q 






v rms 




MM 








10 


2 


0.2 


1000 


0.85 


0.58 


2.35 


8.5 


1.91 X 


10 s 


1.79 x 10 1 


3.66 x 10~ 2 


0.00 


2.07 x 


10 6 


0.86 


10 


2 


0.3 


10 


0.80 


0.53 


2.50 


8.0 


1.07 x 


10 5 


1.37 x 10 1 


3.83 x 10° 


0.46 


2.72 x 


10 6 


0.83 


10 


2 


0.3 


100 


0.88 


0.60 


2.27 


8.8 


2.66 x 


10 5 


1.89 x 10 1 


1.23 x 10~' 


0.00 


1.75 x 


10° 


0.87 


10 


3 


0.2 


10 


0.75 


0.52 


3.98 


7.5 


6.49 x 


10 4 


2.01 x 10 1 


1.77 x 10 1 


0.66 


7.61 x 


10'' 


0.63 


10 


3 


0.2 


100 


1.00 


0.78 


2.99 


10.0 


1.03 x 


10 6 


5.80 x 10 1 


4.43 x 10~ 2 


0.00 


4.41 x 


10" 


0.68 


10 


3 


0.3 


10 


0.83 


0.58 


3.62 


8.3 


1.49 x 


10 5 


2.14 x 10 1 


1.39 x 10 1 


0.57 


6.95 x 


10" 


0.70 


10 


4 


0.1 


10 


0.77 


0.55 


5.18 


7.7 


7.89 x 


10 4 


2.89 x 10 1 


3.37 x 10 1 


0.73 


1.17 x 


10 7 


0.53 


10 


4 


0.1 


100 


0.78 


0.55 


5.15 


7.8 


8.36 x 


10 4 


2.90 x 10 1 


3.34 x 10 1 


0.79 


1.17 x 


1 111 

10 7 


0.54 


10 


4 


0.1 


1000 


0.78 


0.55 


5.15 


7.8 


8.38 x 


10 4 


2.90 x 10 1 


3.33 x 10 1 


0.79 


1.17 x 


10' 


0.54 


10 


4 


0.2 


10 


0.82 


0.59 


4.90 


8.2 


1.31 x 


10 


3.11 x 10 1 


3.14 x 10 1 


0.59 


1.14 x 


10' 


0.58 


10 


4 


0.2 


100 


1.11 


0.92 


3.61 


11.1 


3.26 x 


10° 


1.13 x 10 2 


7.14 x 10~ 3 


0.00 


7.34 x 


10" 


0.61 


10 


4 


0.3 


10 


0.88 


0.63 


4.52 


8.8 


2.80 x 


ur 


3.55 x 10 1 


3.27 x 10 1 


0.53 


1.09 x 


1 r\l 

\{y 


0.64 


12 


2 


0.2 


10 


0.80 


0.51 


2.49 


9.6 


7.47 x 


1 r\A 

l(r 


1.05 x 10 1 


4.58 x 10° 


0.57 


2.80 x 


10° 


0.85 


12 


2 


0.2 


100 


0.90 


0.60 


2.23 


10.7 


2.55 x 


l(r 


1.58 x 10 1 


3.23 x 10~' 


0.00 


1.70 x 


10° 


0.91 


12 


2 


0.3 


10 


0.90 


0.61 


2.22 


10.8 


2.79 x 


l(r 


1.59 x 10 1 


6.30 x 10~' 


0.00 


1.63 x 


I0 ft 


0.91 


12 


2 


0.3 


100 


0.92 


0.62 


2.18 


11.0 


3.34 x 


l(r 


1.68 x 10 1 


4.51 x 10~ 2 


0.00 


1.46 x 


i(r 


0.93 


12 


3 


0.1 


10 


0.82 


0.54 


3.68 


9.8 


8.91 x 


1 /\A 

10 4 


1.79 x 10 1 


1.57 x 10 1 


0.48 


7.11 x 


10 6 


0.67 


12 


3 


0.1 


100 


0.82 


0.54 


3.67 


9.8 


9.03 x 


> /-.A 

10 4 


1.79 x 10 1 


1.56 x 10 1 


0.52 


7.12 x 


10" 


0.67 


12 


3 


0.1 


1000 


0.82 


0.54 


3.67 


9.8 


9.04 x 


1 /\A 

10 4 


1.79 x 10 1 


1.56 x 10 1 


0.54 


7.13 x 


10" 


0.67 


12 


3 


0.2 


10 


0.84 


0.57 


3.57 


10.1 


1.24 x 


10 5 


1.84 x 10 1 


1.36 x 10 1 


0.69 


6.95 x 


10" 


0.67 


12 


3 


0.2 


100 


0.87 


0.60 


3.46 


10.4 


1.78 x 


i(r 


2.00 x 10 1 


1.18 x 10 1 


0.67 


6.66 x 


I0 b 


0.70 


12 


3 


0.2 


1000 


0.89 


0.62 


3.38 


10.6 


2.30 x 


Mr 


2.14 x 10 1 


1.08 x 10 1 


0.76 


6.48 x 


10° 


0.72 


12 


4 


0.1 


10 


0.86 


0.60 


4.65 


10.3 


1.61 x 


Mr 


2.56 x 10 1 


2.64 x 10 1 


0.66 


1.11 X 


\{y 


0.57 


12 


4 


0.1 


100 


0.86 


0.60 


4.63 


10.4 


1.67 x 


Mr 


2.56 x 10 1 


2.61 x 10 1 


0.78 


1.11 X 




0.57 


12 


4 


0.1 


1000 


0.87 


0.60 


4.61 


10.4 


1.78 x 


Mr 


2.56 x 10 1 


2.62 x 10 1 


0.76 


1.11 X 




0.57 


12 


4 


0.2 


10 


0.88 


0.63 


4.53 


10.6 


2.15 x 


10 


2.61 x 10 1 


2.39 x 10 1 


0.63 


1.09 x 


10' 


0.57 


12 


4 


0.2 


100 


0.91 


0.66 


4.38 


11.0 


3.20 x 


10 


2.75 x 10 1 


2.14 x 10 1 


0.65 


1.06 x 


1() 7 


0.62 


12 


4 


0.2 


1000 


1.13 


0.93 


3.53 


13.6 


5.50 x 


10" 


1.33 x 10 2 


1.59 x 10~ 4 


0.00 


7.12 x 


, /id 

10" 


0.61 


12 


4 


0.3 


10 


0.93 


0.67 


4.31 


11.1 


3.96 x 


10 


2.72 x 10 1 


2.17 x 10 1 


0.56 


1.05 x 


10' 


0.66 


14 


2 


0.1 


10 


0.86 


0.55 


2.34 


12.0 


1.13 x 


ur 


9.86 x 10° 


3.49 x 10° 


0.35 


2.23 x 


10° 


0.90 


14 


2 


0.1 


100 


0.86 


0.55 


2.34 


12.0 


1.13 x 


MT 


9.85 x 10° 


3.49 x 10° 


0.36 


2.23 x 




0.90 


14 


2 


0.1 


1000 


0.86 


0.55 


2.34 


12.0 


1.13 x 


Mr 


9.85 x 10° 


3.49 x 10° 


0.36 


2.23 x 


I0 b 


0.90 


1 A 
14 


z 


U.z 


1 n 


u.y 1 


n ^0 
u.jy 


o o 1 
z.Zl 


Iz. / 


Z.4Z X 


J u 


l.Zj X 1U 


L.dL X 1U 


U.J / 


1. 10 X 


1fl6 
LU 




14 


2 


0.2 


100 


0.89 


0.62 


2.24 


12.5 


2.02 x 


10 5 


1.20 x 10 1 


1.17 x 10~ 3 


0.00 


1.36 x 


10 6 


0.85 


14 


2 


0.3 


10 


0.90 


0.62 


2.22 


12.6 


2.26 x 


10 5 


1.24 x 10 1 


7.89 x 10~ 3 


0.00 


1.28 x 


10 6 


0.84 


14 


3 


0.1 


10 


0.89 


0.60 


3.36 


12.5 


2.03 x 


10 5 


1.79 x 10 1 


1.21 x 10 1 


0.51 


6.46 x 


10 6 


0.71 


14 


3 


0.1 


100 


0.89 


0.60 


3.35 


12.5 


2.04 x 


10 s 


1.79 x 10 1 


1.21 x 10 1 


0.54 


6.46 x 


10 6 


0.71 


14 


3 


0.1 


1000 


0.89 


0.60 


3.35 


12.5 


2.05 x 


10 s 


1.80 x 10 1 


1.20 x 10 1 


0.53 


6.47 x 


10 6 


0.71 


14 


3 


0.2 


10 


0.91 


0.62 


3.28 


12.8 


2.72 x 


10 5 


1.91 x 10 1 


1.07 x 10 1 


0.61 


6.29 x 


10" 


0.70 


14 


3 


0.3 


100 


1.08 


0.84 


2.79 


15.1 


3.13 x 


10 6 


7.22 x 10 1 


3.55 x 10" 4 


0.00 


3.59 x 


10" 


0.70 


14 


4 


0.1 


10 


0.94 


0.65 


4.27 


13.1 


3.91 x 


10 5 


2.49 x 10 1 


2.17 x 10 1 


0.63 


1.05 x 


10 7 


0.60 


14 


4 


0.1 


100 


0.94 


0.65 


4.26 


13.2 


4.05 x 


10 5 


2.49 x 10 1 


2.29 x 10 1 


0.48 


1.05 x 


10 7 


0.60 


14 


4 


0.1 


1000 


0.98 


0.67 


4.10 


13.7 


6.98 x 


10 s 


3.11 x 10 1 


2.52 x 10 1 


0.58 


1.02 x 


10 7 


0.61 


14 


4 


0.2 


10 


0.95 


0.67 


4.20 


13.3 


4.90 x 


10 5 


2.55 x 10 1 


2.02 x 10 1 


0.66 


1.04 x 


10 7 


0.59 


14 


4 


0.2 


100 


1.04 


0.74 


3.86 


14.5 


1.74 x 


10 6 


4.71 x 10 1 


2.13 x 10 1 


0.78 


9.60 x 


10" 


0.58 


14 


4 


0.2 


1000 


1.11 


0.79 


3.59 


15.6 


5.46 x 


10 6 


8.42 x 10 1 


3.14 x 10 1 


0.84 


9.15 x 


10" 


0.58 



